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using three dimensional design techniques and the 
Advanced Surface Design software on the Computers ion 
Designer V-X Interactive Graphics system, the aircraft 
configuration was created. The canard, tail, vertical 
tail, and main wing were created on the system using 
'Wing Generator-, a computervision based program 
introduced in Appendix A. 2. These plots can be seen in 
Figures A. 1.1 - a. 1.3. The individual components of the 
plane were created separately and were later 
individually imported to the master database. An 
isometric view of the final configuration can be seen 

in Figure A. 1.4. 
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INTRODUCTION 


The creation of a complex entity such as an 
airplane wing can be a complex and time consuming task 
on a CAD system. This program allows the designer a 
simpler, faster, automated means by which to experiment 
with wing designs and create wing models. These models 
may then be used for creating plots, defining mass 
properties, and doing aerodynamic and structural 
analyses as necessary. 

In general, a wing's airfoil or cross section is a 
function of the spanwise location. For this reason, an 
input file may be created which defines up to ten 
airfoils to be used along the span of the wing. If one 
does not wish to create airfoils (ie. a simple wing 
geometry), a default file containing symmetrical 
airfoils may be utilized instead. At spanwise 
locations where an airfoil section is not specified, a 
linear interpolation provides a smooth transition of 
the wing surface. 

The program format requires that the user input 
some basic data on the dimensions of the wing; it then 
responds with a CAD model of the wing, and provides 
plots of all views of the wing. 
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PROGRAM THEORY 


A basic design assumption for wing Generator 1.0 is 
a constant taper ratio, defined by the lengths of the 
tip and root chord airfoils. Once these parameters are 
defined, the chord length of any intermediate airfoil 
is known; it is a function of the span location only. 

As airfoils in the input file are of length unity, they 
may simply be scaled by the taper ratio at the required 

span location. 

After being scaled to the correct size, the 
coordinates of each airfoil are adjusted to reflect the 
angle of twist defined by the user. At this stage, the 
cartesian coordinates of the airfoil sections are stored 
in a database; they are then connected together using 
the CAD BSPLINE entity and copied into model space to 
form a complete wing. Finally the program interacts with 
the user to detail the wing and plot different views. 

WING DESIGN ILLUSTRATION 

in order to demonstrate the graphic capabilities of 

the program, a wing was designed with an elliptic span 

load distribution. It was to be a forward swept wing, 

with the following characteristics: 

: 50 ft. 

: 750 sq. ft. 
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Wing Span 
Wing Area 



Root Chord : 20 ft. 

Tip Chord : 10 ft. 

Sweep Angle : -30 degrees 
Lift Coeff. : 1.00 

No. Airfoils: 6 . ^ j f *. 

Knowing the number of airfoils and total lif 

coefficient for the wing, an expression was derived for 

the local spanwise lift coefficient, (C^). For any 

elliptically loaded wing, the load at any spanwise 

position is given by equation (A. 2.1). (See Fig. A. 2.1) 


Elliptic Load Distribution 



w(y)=w 0 yi-(y/<*>/2j7 T 1 

However, the load at the root is given by eq. (A. 2. 2). 
W 0 =C L ^V 2 S(4/*b) (A. 2. 2 
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and the total lift on 


the wing by (A. 2. 3) 


c(y)-c 1 (y)dy*sfV J 

0 

integrating eg. (A. 2. 3) yields 


(A. 2. 3) 


l= c l*^ v ’ s 


(A. 2. 4) 


and 


b/2 


C L = (2/s) 


c(y) •c 1 (y)dy 


(A. 2. 5) 


Substituting 


equation (A. 2. 2) into equation (A. 2.1) yields 


w(y)=c L -^pv I s-(4/irb)/i-(y/(b/2)) / 


(A. 2. 6) 


Substituting equation (A. 2. 5) into equation (A. 2. 6), 


Dividing both Eq. (A.2.6) and (A.2.7) by the dynamic pressure 
and solving tor the local spanwise lift coefticient, 

c,(y)-c,.».(4/«.).A35aBZain (A. 2. 8) 

1 L C(y) 

Thus for the wing described above, the local lift 
coefficient may be given as eq. (A.2.8), provided c(y) 

is given by: 
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(A. 2.9) 


c(y)=c r ~(y/(b/2) ) [c r -c t ] 

where C r is the root chord 
and C t is the tip chord. 

Substituting the values for this wing into eg. (A 
the following relation is found: 

Cl (y)=l9.10-^=SZ2lII (A 

1 ( 20 - 0 . 4y) 

The data for the six airfoils as derived from eq. 
(A. 2. 10) are given below: 


Table A. 2.1 

Span Locat ion (v) 
0 
5 

10 

15 

20 

25 


Tift, coe fficient (Cll 
0.955 
1.040 
1.094 
1.091 
0.955 
0.000 


. 8 ) , 


2 . 10 ) 
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Uft Coefficient 


Figure A. 2. 2 


Lift Coefficient Vs. Span 



Wing generator 1.0 requires actual airfoil 
coordinates as input. Thus an intermediate step is 
required, that of creating airfoil section coordinates 
from a target lift coefficient. The program 
'Design. for' (Reference 13) was utilized for this task, 
the task of creating some simple airfoils of varying 

geometry along the wing span. 

'Design. for' uses target pressure distributions, given 

as input, to determine airfoil section coordinates. 
Examples of input pressure distributions and format of 
output are shown at the end of this appendix. 
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Target pressure distributions were prescribed which 
would yield a value in Table (A.2.2) for that airfoil. 
(The output from 'Design. for' for these airfoils may be 

found at the end of this appendix) . 

Finally, the airfoil coordinates from 'Design. for’ 

were merged and formatted as input to the wing 
Generator. (The method for transferring this file from 
the PC to the CAD system is described in this appendix 
(see 'File Transfers'). The program was run for the 
parameters of the forward-swept wing and airfoils 
described previously; plots are included at the end of 

this appendix. 
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USING WING GENERATOR 1.0 


Wing generator 1.0 is a fully documented computer 
program, written with the aircraft designer in mind. 

This step by step guide should aid the first time user 
to run the program successfully. 

computer: Please enter the date. 
user: Type the date. 

mmnuter: Please enter the wing span. 

„ ser: Type the length of the wing from tip to tip. 
renter: Please enter the length of the root chord. 
„ sgr: Type the length of the airfoil located at the 

fuselage . 

renter: Please enter the tip chord length. 

„ ser; Type the length of the airfoil furthest from the 

fuselage . 

reenter: Please enter the sweep angle in degrees. 
user; Ty pe the sweep angle using the convention that a 
positive angle indicates a rearward swept wing. 
reenter: What is the airfoil section filename? 

(Default = ( P . THJ . FOIL) 

„ ser; Hit the enter key to choose the default airfoil 
section, or type in another previosly created 
data file. The default file is made up of simple 
symmetrical airfoils. The procedure to make up 
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new data files will be covered in the section 
titled File Transfers. 

computer: There are: X airfoils contained in the file. 

Airfoil section (1) has been uploaded. 

At how many span stations will (1) appear? 
user: The computer has now read in all "X" of the 

airfoils. It now runs through each airfoil from 
root to tip and asks you how many times each of 
these airfoils will be used to define the wing 
geometry. Type the appropriate number to 
continue . 

computer: Please indicate the # (1) span position for 
airfoil 1. 

Span position =? 

user : Type in the the location of the airfoil. The 

location is defined as the perpendicular distance 
between the airfoil under consideration and the 
root airfoil. 

router: Please enter the twist angle in degrees for 
same . 

Twist angle =? 

user: Type the number of degrees that the airfoil under 
consideration is twisted. This angle is equal to 
the negative of the angle of attack. 

NOTE : This routine of defining the position and twist 
angle for each airfoil is repeated until all of 
the airfoils in the data file are defined. 
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^n„ter: Do you want to see the airfoil coordinates? 
user: Hit enter to default a "no" answer or type "y" to 
see the true space coordinates printed. 
ffimmitfir : DO you want to see the points plotted? 

„ ser; Hit enter to default a "no" answer or type "y" to 
create points at each of the coordinates along 

the airfoil. 

rnwniiter: DO you want to surface the wing? 

„ ser; Hit enter to default a -no” answer or type ”y” to 
create a surface on the wing geometry. 
rnmnnter: Activate Part <var> 

aserj. Type the name you wish the part to have. Hit the 
"Ctrl" and "x" keys simultaneously. 
f-nmnuter : Dynamic view speed 5: view 

„ gar: use the pen to digitize a view, then use the ICU 
to adjust to the location that looks best. Hit 

Ctrl x. 

^„n.,ter: Insert label "Leading edge” : Draw/Model entity 

„ g „ r; select the modifier "near" with the pen and 

digitize the leading edge line and the location 
for the label to go. Hit ctrl x. 

^n.nnter: Insert label "Trailing edge": Draw/Model 

entity 

user: Repeat last step for trailing edge line. 
^■nuter: Insert label "Quarter chord" : Draw/Model 

entity 

user: Repeat again for quarter chord line. 
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r.nmnuter : Insert Bspline Tag = HI:Model location 
user: Select the modifier "near" with the pen and then 
digitize the upper points on the airfoils as 


BEFORE COUUANC 



FRONT VIEW OF WING 

AFTER COUUANL 




WP! CAD LABORATORY 

TITLE. FIGURE i 

NO. I 

SCALE.0.0 DATE. 

DRAWN BY: tj a NPf 

SHEET, i 


rnmnuter : Insert label "Upper Boundary" : Draw/Model 
entity 

user: Repeat previous procedure for inserting labels 
for the upper boundary line. 
rnmnuter: Insert Bspline Tag = LOW:Model location 
user: Select the modifier "near" with the pen and then 
digitize the lower points on the airfoils as 


shown in Fig. 

2. Hit Ctrl x. 



BEFORE COMUAHO 


Cl ' 

1 "i i 1 

K D3 ^ 

<=> 


r RON 

T VIEW OF WING 



after couwanc 



1 


WT! CAD LABORATORY 

TITLE. FiCORE 2 

NO. l 

SCALE.C 0 DATE. W/50 

DRAWN BY: thj a m 

SHEET, i 






insert Ntext '"'wing section»THJ t NPF»scale"l 
"date" 1 " Draw ent 

user: Use the pen to digitize any point on the title 
block. Hit Ctrl x. This procedure will 
automatically configure the title block with the 
appropriate scale and date. 


In order to have multi-color plots, entities were 
created on different layers in the drawing. The 
entities and their respective layers are as follows: 


Airfoils Layer 0 

Title block Layer 0 

Leading Edge Layer 1 

Trailing Edge Layer 2 

Quarter Chord Layer 3 

Upper Bound Layer 4 

Lower Bound Layer 5 

Rotated Airfoils Layer 6 
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CONCLUSION 


It has been shown that the time required to create 
and modify airplane wings using the Coroputervision 
Interactive Graphic System has been reduced 
significantly. The program Wing Generator 1.0 should 
prove to be an invaluable tool to the airplane 
designer, as it quickly and efficiently creates wing 
geometries which may be used as models for structural, 
aerodynamic, and mass properties analysis. 


A. 2. 14 


forward swept wing statistics 


Wing Span 

Root Chord Length 
Tip Chord Length. 
Sweep Angle 


100 Ft. 
.20 Ft. 
.10 Ft. 
-30 Deg. 


g pan station 
1 
2 

3 

4 

5 


ppan Loc ationf Ft Y 
0 

12.5 

25.0 

37 . 5 

50.0 


Twist An ql e ( Deg ) 
-5 

-2.5 

0 

2.5 

5 
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INPUT FILE FORMAT 

The airfoil data file may contain up to ten 
airfoils, each being made of fifty distinct points (25 
upper, 25 lower). Each airfoil in the input file must 
have a chord length of unity, and twist angle of zero 
degrees, with the leading edge point coinciding with 
the origin (0,0) and the trailing edge point coinciding 
with (1,0). Wing Generator is formatted to accept 
output from the program 'Design. for' (Reference 13). 

The first entry in the file should be an index which 
represents the number of airfoils contained in the 
file. A sample input file is shown in this appendix. 


Input File Format: . . , 

- number of airfoils in file 

XXXXXXX YYYYYYY - coordinates of first point 

— rest of coordinates, 1st foil 

XXXXXXX YYYYYYY - coordinates of last point 
XXXXXXX YYYYYYY - 1st point, 2nd foil coordinates 
- etc. 
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Wing 


Generator Sample Input File 


1 

.000000 
. 004278 
.017037 
.038060 
. 066987 
.103323 
.146447 
.195619 
. 250000 
. 308658 
.370590 
.434737 
. 500000 
.565263 
.629410 
.691342 
.750000 
.804381 
.853553 
.896677 
.933013 
.961940 
.982963 
.995722 
1.000000 
.995722 
.982963 
.961940 
.933013 
.896677 
.853553 
.804381 
.750000 
.691342 
.629410 
.565263 
.500000 
.434737 
. 370590 
.308658 
. 250000 
.195619 
. 146447 
.103323 
.066987 
.038060 
.017037 
.004278 
.000000 


.000000 
.017986 
.041356 
.065242 
.089596 
.112707 
.134173 
.152904 
.168457 
.180038 
.187144 
. 189127 
. 184654 
.172402 
.154952 
.134813 
.113580 
.092567 
.072772 
.054881 
.039328 
.026240 
.015517 
.006885 
.000000 
-.005329 
-.009038 
-.011078 
-.011761 
-.011726 
-.011856 
-.012942 
-.015570 
-.019980 
-.026221 
-.033723 
-.041094 
— . 046654 
-.050104 
-.051724 
-.051922 
-.050678 
-.048149 
- . 044114 
-.038737 
-.031603 
-.022988 
-.012153 
.000000 
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file transfers 


Generally data files defining the airfoils will be 
generated outside the Computervision system. The 
program being used in this project is 'Design. for' 
(Reference 13). In these cases, the output file of the 
source program must be in the form acceptable to Wing 
Generator 1.0, described previously. Once the file is 
in the correct format and saved on a floppy disk, it is 
ready to be transferred. The following procedure 
should aid the user in completing such transfers 
successfully • 

NOTE ; The user should type what is UNDERLINED - 

1 . on the patch panel in the CGP room (back room 
of the lab) , connect a jumper cable between 

« CV9600 BAUDPCU" and "IBM XT#1" (see proctor 
for help) . 

2. Boot up the IBM XT. Type the correct date and 
time. Do not put a floppy diskette in the 
drive when booting up. Make sure the T 
switch" the IBM is on "Patch panel". At the 
main menu, choose the option to get to DOS 

level . 

3. The program Kermit is used to transfer (ASCII) 
data from the IBM XT to the Computervision 
CADDS 4X system or vice versa. Enter the 
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Kerniit program by typing 
C:\>KEBMIT <CR> 

IBM-PC Kermit-MS V2.26 
Type ? for help 

Kermit-MS> (This is IBM Kermit Prompt) 
To set the necessary default parameters, type 
Kermit-MS>TAXF, gFT.CAg <CR> 
Kermit-MS>£EL_BAUD_9M0 <CR> 

Place the floppy disk you want the data file 
transferred from in drive (A). Specify that 
the file be transferred from the floppy disk 
drive by typing 

Kermit-MS>SE2LJ2EFA1ILT_A1 <CR> 

TO connect the IBM to the CV system, type 

Kermit-MS>CONHECT <CR> 

On the CV system, type 
1 > KERMIT <CR> 

Kermit server(This is CV Kermit Prompt) 
At this time press the EXE1 and i keys 
simultaneously to log onto the CV system. Type 
your initials at the prompt 
TYPE NAME, NUMBER 
XXX <CR> (Your initials) 

** TASK # INITIATED ** 

At this time return back to Kermit on the IBM 
by pressing the CT RL and 1 key simultaneosly 
and then the C key (no carriage return). The 



10 . 


11 


12 


13. 


14 


15 


IBM Kermit prompt "Kermit-MS>" will appear. 

The command to begin the transfer procedure is 
rrrmi t-MS> £EM& XXXXXXXX-XXi i <CR> 
once the file has been transferred, type 
Kermit-MS>CONNECX <CR> 

This will return the user to the CV Kermit 
prompt -Kermit server". To exit Kermit on the 
CV system, press the CTRL key and then the C 

Key twice: 

Kermit server CTRL CC 

Once at the CV system level prompt, log out by 
typing 

1> LQG <CR> 

once again return to the IBM Kermit prompt by 
typing 

Ctrl 1C (no carriage return) 

Exit Kermit on the IBM by typing 
Kermit-MS>£XIT <CR> 

C:\> 

nt transfer is complete. 


The original filename (first level) should not be 
more than eight oharaters in length and the catalog 
level (second level) not more than three characters. 
This file must now be renamed on the CV system to a 
multilevel structure so that it may be manipulated by 
daily save procedures. To rename the file to this 
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multilevel structure, the COPYTEXT command may be used 
at system level as follows: 

1 > COPYTEXT XXXXXXXX.XXX C . HI , xxx xXXXX • XX % 

Where "III" represents the user's initials. 
Once these steps are completed, the user can run 
Wing Generator. When the program asks for the data 
filename, enter the new CV system name. 
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Sample Input File to 'Design. for' Program 


NT = 

25 

NCP = 

22 

01 

0.00 

02 

0.010 

03 

0.025 

04 

0.050 

05 

0.075 

06 

0.10 

07 

0.40 

08 

0.45 

09 

0.50 

10 

0.55 

11 

0.60 

12 

0.65 

13 

0.70 

14 

0.75 

15 

0.80 

16 

0.85 

17 

0 . 90 

18 

0.925 

19 

0.950 

20 

0.975 

21 

0.990 

22 

1.000 


- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.000 

- 1.000 

+ 0.012 

- 0.960 

+ 0.050 

- 0.765 

+ 0.110 

- 0.567 

+ 0.183 

- 0.408 

+ 0.238 

- 0.290 

+ 0.270 

- 0.205 

+ 0.293 

- 0.137 

+ 0.300 

- 0.085 

+ 0.285 

- 0.043 

+ 0.235 

- 0.027 

+ 0.190 

- 0.015 

+ 0.142 

- 0.005 

+ 0.068 

- 0.001 

+ 0.020 

+ 0.000 

+ 0.000 Sample 
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Output from 'Design. for' Program 


DESIGN. FOR. .M. S .GARELICK. .06-09-89 
MILNE-THOMSON THIN AIRFOIL DESIGN PROCEDURE 
THICKNESS. .CAMBER. .COORDINATES 
NT=2 5 NCP=2 2 


INPUT CP DISTRIBUTION 
XCP= .500000 CPU= -1.193000 CPL= 

X CP= .490000 CPU= -1.193000 CPL= 

XCP= .475000 CPU= -1.193000 CPL= 

X CP= .450000 CPU= -1.193000 CPL= 

XCP= .425000 CPU= -1.193000 CPL= 

X CP= .400000 CPU= -1.193000 CPL= 
XCP= .100000 CPU= -1.193000 CPL= 
XCP= .050000 CPU= -1.193000 CPL= 
XCP= .000000 CPU= -1.193000 CPL= 
XCP= -.050000 CPU= -1.073000 CPL= 
XCP= -.100000 CPU= -.954000 CPL= 
XCP= -.150000 CPU= -.835000 CPL= 
XCP= -.200000 CPU= -.716000 CPL= 
XCP= -.250000 CPU= -.596000 CPL= 
XCP= -.300000 CPU= -.477000 CPL= 
XCP= -.350000 CPU= -.358000 CPL= 
XCP= -.400000 CPU= -.239000 CPL= 
XCP= -.425000 CPU= -.179000 CPL= 
XCP= -.450000 CPU= -.119000 CPL= 
XCP= —.475000 CPU= —.060000 CPL= 
XCP= -.490000 CPU= -.024000 CPL= 
XCP= -.500000 CPU= .000000 CPL= 


.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.000000 

.012000 

.050000 

.110000 

.183000 

.238000 

.270000 

.293000 

.300000 

.285000 

.235000 

.190000 

.142000 

.068000 

.020000 

.000000 


CL= 1.000210 
1=25 XA= .500000 PC= 
1=24 XA= .495722 PC= 
1=23 XA= .482963 PC= 
1=22 XA= .461940 PC= 
1=21 XA= .433013 PC= 
1=20 XA= .396677 PC= 
1=19 XA= .353553 PC= 
1=18 XA= .304381 PC= 
1=17 XA= .250000 PC= 
1=16 XA= .191342 PC= 
1=15 XA= .129410 PC= 
1=14 XA= .065263 PC= 
1=13 XA= .000000 PC= 
1=12 XA= -.065263 PC= 
1=11 XA= -.129410 PC= 
1=10 XA= -.191342 PC= 
1= 9 XA= -.250000 PC= 
1= 8 XA= -.304381 PC= 
1= 7 XA= -.353553 PC= 
1= 6 XA= -.396677 PC= 
1= 5 XA= -.433013 PC= 
1= 4 XA= -.461940 PC= 
1= 3 XA= -.482963 PC= 


-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.596500 PT= -.596500 
-.600668 PT= -.592332 
-.621500 PT= -.571500 
-.584479 PT= -.452195 
-.549678 PT= -.334327 
-.500533 PT= -.236074 
-.444500 PT= -.151500 
-.382630 PT= -.083944 
-.315495 PT= -.034048 
-.242616 PT= -.004293 
-.167193 PT= .007423 
-.098740 PT= .007918 
-.041704 PT= .000815 
A. 2. 23 


1= 2 XA= -.495722 PC= -.009411 PT- - 
1= 1 XA= -.500000 PC= .000000 PT- 
ALFADEG= 2.033764 

CAMBERLINE SLOPES AND COORDINATES 

.500000 CSLOPE= -.776462 ZC 
.495722 CSLOPE= 

.482963 CSLOPE= 

.461940 CSLOPE= 

.433013 CSLOPE= 

.396677 CSLOPE= 

.353553 CSLOPE= 

.304381 CSLOPE= 

.250000 CSLOPE= 

.191342 CSLOPE= 

.129410 CSLOPE= 

.065263 CSLOPE= 

.000000 CSLOPE= 

-.065263 CSLOPE= 

-.129410 CSLOPE= 

-.191342 CSLOPE= 

•.250000 CSLOPE= 

-.304381 CSLOPE= 

•.353553 CSLOPE= 

-.396677 CSLOPE= 

-.433013 CSLOPE= 

-.461940 CSLOPE= 

-.482963 CSLOPE= 

-.495722 CSLOPE= 

-.500000 CSLOPE= 

,268364 TAU= 


,000856 

.000000 


1=25 XA= 
1=24 XA= 
1=23 XA= 
1=22 XA= 
1=21 XA= 
1=20 XA= 
1=19 XA= 
1=18 XA= 
1=17 XA= 
1=16 XA= 
1=15 XA= 
1=14 XA= 
1=13 XA= 
1=12 XA= 
1=11 XA= 
1=10 XA= 

1= 9 XA= - 
1= 8 XA= 

1= 7 XA= 

6 
5 
4 
3 
2 
1 


1= 

1 = 

1 = 

1 = 

1 = 

1 = 


XA= 

XA= 

XA= 

XA= 

XA= 

XA= 


ALAMBDA= 


-.587146 ZC= 
-.395264 ZC= 
-.331109 ZC= 
-.264196 ZC= 
-.223868 ZC= 
-.180330 ZC= 
-.149182 ZC= 
-.113915 ZC= 
-.086898 ZC= 
-.054005 ZC= 
-.030686 ZC= 
.014020 ZC= 
.060764 ZC= 
.094334 ZC= 
.130074 ZC= 
.156723 ZC= 
.181352 ZC= 
.199114 ZC= 
.212741 ZC= 
.216285 ZC= 
.212537 ZC= 
.200456 ZC= 
.185470 ZC= 
.178115 ZC= 
.094476 


FOURIER COEFFICIENTS 
1= 1 AN= . 318473D+00 BN= .173889D+00 
1= 2 AN=- . 920674D-01 BN= .528346D-01 
1= 3 AN= . 390524D-01 BN=- • 255276D-01 
1= 4 AN=- • 598885D-01 BN=- . 269270D-01 
1= 5 AN= . 342235D-01 BN= .557217D-03 
1= 6 AN=- . 318984D-01 BN= .840308D-02 
1= 7 AN= . 208790D-01 BN=- . 110097D-02 
1= 8 AN=- . 238717D-01 BN=- . 453781D-02 
1= 9 AN= . 191965D-01 BN= .108491D-02 
1=10 AN=- . 161761D-01 BN= .344226D-02 
1=11 AN= . 127289D-01 BN=- • 183084D-03 
1=12 AN=- . 128800D-01 BN=- . 207360D-02 
1=13 AN= . 117038D-01 BN= .285786D-03 
1=14 AN=- . 964384D-02 BN= .176353D-02 
1=15 AN= . 723877D-02 BN= .802936D-04 
1=16 AN=- . 722481D-02 BN=- . 117083D-02 
j = 17 AN= .701074D-02 BN= .689152D-04 
1=18 AN=- . 518064D-02 BN= .112414D-02 
1=19 AN= . 338631D-02 BN= .977638D-04 
1=20 AN=- . 321564D-02 BN=- . 899677D-03 
1=21 AN= . 339589D-02 BN=- . 958598D-04 
1=22 AN=- • 163063D-02 BN= .744816D-03 
1=23 AN=- . 317677D-04 BN=- . 307924D-04 

A. 2. 24 


.000000 

.002916 

.009184 

.016819 

.025430 

.034297 

.043012 

.051113 

.058267 

.064157 

.068520 

.071236 

.071780 

.069340 

.064365 

.057416 

.049005 

.039812 

.030458 

.021578 

.013783 

.007581 

.003240 

.000778 

.000000 



1=24 AN= . 142627D-08 BN=- . 817194D-03 
1=25 AN= . 317667D-04 BN=- . 307941D-04 
THICKNESS SLOPES AND COORDINATES 


1=25 XA= 

.500000 TSLOPE= 

-.134182 

ZT= 

1=24 XA= 

.495722 TSLOPE= 

-.132822 

ZT= 

1=23 XA= 

.482963 TSLOPE= 

-.128734 

ZT= 

1=22 XA= 

.461940 TSLOPE= 

-.122044 

ZT= 

1=21 XA= 

.433013 TSLOPE= 

-.112648 

ZT= 

1=20 XA= 

.396677 TSLOPE= 

-.100999 

ZT= 

1=19 XA= 

.353553 TSLOPE= 

-.086625 

ZT= 

1=18 XA= 

.304381 TSLOPE= 

-.070540 

ZT= 

1=17 XA= 

.250000 TSLOPE= 

-.051456 

ZT= 

1=16 XA= 

.191342 TSLOPE= 

-.031289 

ZT= 

1=15 XA= 

.129410 TSLOPE= 

-.006589 

ZT= 

1=14 XA= 

.065263 TSLOPE= 

.019587 

ZT= 

1=13 XA= 

.000000 TSLOPE= 

.057739 

ZT= 

1=12 XA= 

-.065263 TSLOPE= 

.089432 

ZT= 

1=11 XA= 

-.129410 TSLOPE= 

.097954 

ZT= 

1=10 XA= 

-.191342 TSLOPE= 

.099756 

ZT= 

1= 9 XA= 

-.250000 TSLOPE= 

.094058 

ZT= 

1= 8 XA= 

-.304381 TSLOPE= 

.085446 

ZT= 

1= 7 XA= 

-.353553 TSLOPE= 

.074532 

ZT= 

1= 6 XA= 

-.396677 TSLOPE= 

.064434 

ZT= 

1= 5 XA= 

-.433013 TSLOPE= 

.056483 

ZT= 

1= 4 XA= 

-.461940 TSLOPE= 

.051379 

ZT= 

1= 3 XA= 

-.482963 TSLOPE= 

.048758 

ZT= 

1= 2 XA= 

-.495722 TSLOPE= 

.047660 

ZT= 

T = 1 XA= 

-.500000 TSLOPE= 

.047238 

ZT= 

AIRFOIL COORDINATES. . XUPPER ( I ) =XLOWER ( I ) 

.000000 

.000000 



.004278 

.017986 



.017037 

.041356 



.038060 

.065242 



.066987 

.089596 



.103323 

.112707 



.146447 

.134173 



.195619 

.152904 



.250000 

. 168457 



.308658 

.180038 



. 370590 

.187144 



.434737 

.189127 



.500000 

. 184654 



.565263 

. 172402 



.629410 

.154952 



.691342 

.134813 



.750000 

.113580 




.804381 

.092567 


.853553 

.072772 



.896677 

.054881 



.933013 

.039328 



.961940 

. 026240 



.982963 

.015517 



.995722 

.006885 



1.000000 

.000000 




.000000 
.017475 
.034594 
.051008 
.066368 
.080351 
.092631 
.102918 
.110902 
.116318 
.118797 
.117946 
. 112885 
.103253 
.090988 
.078048 
.065363 
.053615 
.043144 
.034049 
.026135 
.019076 
.012522 
.006211 
.000000 
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.995722 

.982963 

.961940 

.933013 

.896677 

.853553 

.804381 

.750000 

.691342 

.629410 

. 565263 

. 500000 

.434737 

. 370590 

. 308658 

.250000 

.195619 

.146447 

. 103323 

.066987 

. 038060 

.017037 

.004278 

.000000 


-.005329 

-.009038 

-.011078 

-.011761 

-.011726 

-.011856 

-.012942 

-.015570 

-.019980 

-.026221 

-.033723 

-.041094 

-.046654 

-.050104 

-.051724 

-.051922 

-.050678 

-.048149 

-.044114 

-.038737 

-.031603 

-.022988 

-.012153 

.000000 



PROGRAM LISTING 


WING GENERATOR 1.0 

WRITTEN BY 

TOM JUTRAS & 
NOAH FOR DEN 

1 ? 89- 1 990 


<# 

<0 
<# 

<0 

DIM X (50) » Y (50) »FOX (2000) »FOY <2000) 

*( £ *”“* 

<0 
<0 
<0 
<0 
<0 
<0 
<0 
<0 
<0 

<0 

MAXCOR = 49 
INC - 0 

SFNAME = "P.THJ. FOIL- 

ACT PART <VAR> „ TP ... , 

ACT DRA DETAIL! FORM L.KF.S.C DRAW 4 

SEL CPI. 3 
REG GRA TAG 
SEL TAG ON 

F’RNT PLEASE ENTER THE WING SPAN IN INCHES 
PR NT 

READ (WING SPAN -? ) B 

f'rnt please enter the length or THE root chord in inches 

READ (ROOT CHORD LENGTH ? ) RL 

PRNT PLEASE ENTER THE TIP CHORD LENGTH IN INCHES 
READ (TIP CHORD LENGTH s’ ) TC 

PRNT PLEASE ENTER THE SWEEP ANGLE IN DEGREES 
PRNT 

READ (SWEEP ANGLE *? )SWEfcr- 
PRNT 

XMAX = ( B/Z )* TAN (SWEEP )+TC 

CHA L VIE scale tscalemname top* name ISO! name front INAME 5 

PRNT WHAT IS THE AIRFOIL SECTION FILENAME (DEFALL 

READ &FNAME 
OPE NR 1 » SFNAME 
READF 1 • &NUMSC 
NUMSFC - SNUMSC (1 >7) 

PRNT THERE ARE • -f NUMSEC) AIRFOILS CONTAINED IN I HE FILE 

COUNT = 1 
REPEAT 


RIGHT 

iFNAMEJ) 
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NODE = 1 
REPEAT 

READF 1* &TXT 
X (NODE ) = 61 XT ( 1 » 3 0> 

Y (NODE) = &TXT < 3 3 » Zif > 

NODE - (NODE + !>■ 

UNTIL (NODE . EO . (MAXCQR + 1 '> '> 

S i&jm w» wwar*-» 

READ (TCvANST) LANS 
ANS = 0. A N S (1*7) 

INDEX = 1 
REPEAT 

.( ft 

■(# GET ALL PARAMETERS HERE 

sb PLEASE JNBJCATE THE f «!«•«» SPAN POSITION FOR AIRFOIL <COUN 

PRINT , „ , „ 

READ (SPAN POSITION 

prn T t please enter the twist angle in degrees for same 

READ (TWIST ANGLE *? ) TWIST 

<#REST OF ROUTINE to INSERT AIRFOILS C-ULS HERE 
<# 

XLE = 7-* TAN (SWEEP) 

XTE= ( (XMAX-RC) / (P/2) )*2+RC 

CLEN-XTE-XLE 

XQUAR-XLE+ ( ,25*CLEN) 

NODE = & 

REPEAT 

NODE = (NODE + 1) 

SPOT = (<10<T * INC) + NUDE^ 

FOX (SPOT) = ( X (NODE) * GLEN) +XLE 
FOY (SPOT) - Y (NODE) * CLEN 
TEMPX * - FOX (SPOT) 

TEMPT - FOY (SPOT) 

if' (TEMPX .LT. XOUAR) ADD = . , TEf!PX - XOUAR) ) ) 

J„AR ?I1?E^X^RM ♦ <TEMPT * TEMPT,. “ *•» 
m (ipm, ■ (WEC'S (COS (ARC,) - XOUAR 
FOY (SPOT) = (RVEC * (SIN (ARt-O ) ) 

UNTIL (NODE . EQ. MAXCUR) 

FOX (SPOT +D - 1 
FOX (SPOT + Z) - TWIST 
FOX. (SPOT + 3) - XQUAR 
INC = (INC + 1) 

INDEX = (INDEX + 1) 

UNTIL (INDEX .EQ. (ANS+D) 

COUNT = (COUNT + D ORIGINAL PAGE IS 

n'f PDca 0 l ALiTY 



UNTIL (COUNT .EQ. (NUMSEC +1)) 

■(# 

wmoHs TwSS? MS? 1 *? <Z> 

<*$£& S£ S8i?iS 11 ^ « cSoRO POSITION AT <Z> 

<# 

REftJ (BOTTOM WANT TO SEE THE AIRFOIL COORDINATES ?> SLST 
IF (&LST -EQ. "N") GOTO JUMP 

PR NT 
PRNT 

mil —AIRFOIL LIST— TRUE SPASI-I COORDINATES. 

PRNT 
POL = 0 
REPEAT 

pfil - t POL + ) ) 

PRN'f {pi'll >» \ FOX (POL ) y > \ F OY (POL) > 

UNTIL (POL .EQ. (INC * 100) > 

# JUMP 

< # 

it ROUTINE TO PLOT D8PL1NES OF AIRFOILS IN MASTER ARRAY 

< # 

<# 

NUN = 0 
ZNUM = 50 
CHECK =48 
REPEAT 

INS ESPL ! '•# 

REPEAT 

NUN = NUM + 1 „ 

XC FOX (NUM) y'ii FOY (NUM) 37TFUX w.NUM) 1 

IJNT I L ( NUM . GT . CHECK ) 

<CR> 

NUM = NUM + 51 

CHECK = CHECK + 100 

ZNUM = ZNUM +100 

UNTIL (ZNUM .G1. (INC* 100)) 

DYN VIEW SPEED 5! <VAR> 

IeAP (b1 N U.U WANT TO SEE THE POINTS PLOTTED ?) SAXE 

IF (&AXE . E Q . " N ' ' ) GOTO LEDGE 

I NIT FEM 

TOM = 0 

REPEAT 

NOAH = 1 

REPE lJ s GPOl :X-!FOX (NOAH 4 (TOM* 100) ) >Y< FOY (NOAH * <T0H«1S0>>> 

Z-f FOX (T DM* 100 + 50)> 

NOAH - NOAH +5 
UNTIL (NOAH .GT. MAX COR) 
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TOM - TOM + 1 
UNTIL (TOM .EG, INC) 

<# 

l# THIS ROUTINE WILL CONNECT THE LEADING EDGES 

<# 

<# 

#LEDGE 
SEL LAY 1 
NUM = 1 

CHECK = JNC*100 

INS BSF'L TAG LEDGE: <# 

X-CFOX (NUM) 3YTFOY (NUM) DZI FOX (NUM + 49)}><# 

NUM = NUM + 3 0$ 

UNTIL (NUM .GT* CHECK) 

< CR > 

INS LAB “LEADING EDGE 11 SNEAK *\VAR/ 

<# 

■<# THIS ROUTINE WILL CONNECT THE TRAILING EDGES 

<# 

<# 

SEL LAY 2 

NUM = 25 

CHECK = INC-* 100 

INS BSF’L TAG TEDGE : <* 

XT FOxIniJM ) mr-OY ( NUM ) D 7A F OX ( NUM + 25) 5 y < # 

NUM -• NUM + 3 00 
UNTIL (NUM . GT. CHECK) 

.( f ft \ 

INS LAB "TRAILING EDGE" '■ NEAR <VAR> 

<# 

<* THIS ROUTINE KILL CONNECT THE QUARTER CHORD POINTS 

<# 

<# 

SEL LAY 3 

NUM = 52 

CHECK = INC* 100 

I N S B S P !.. T A G G U A R C H : # 

REPEAT 

X \ FOX ( NUM ) ) Y 97i F UX ( NUM • o )>■-.£ 

NUm’= NUM +100 
UNTIL (NUM . GT. CHECK) 

INS LAB "QUARTER CHORD": NEAR <VAR> 

ZOOM DP A WIN: <VAR> 

SEL LAY 4 

INS BSPL TAG Hi: NEAR <VAR> 

INS LAB "UPPER BOUNDARY" : NEAR <VAR,> 

SEL LAY 5 

INS esp |- TAG LOW: NEAR \ VAR ? 
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INS LAB "LOWER BOUNDARY NEAR <VAR> 

ZOOM DRA ALL 

<# 

<# THIS ROUTINE WILL DETAIL THE DRAWINGS 

<# 


PRASE ENT :V1E NAME TOP TAG HI LOW 
ERASE ENT :VIE NAME ISO TAG HI LOW GUnRCH 
ERASE ENT :V3E NAME RIGHT TAG’ QUARL-H 
ECHO TAG OFF 
ft EL MOD DRA 

or An (PI EASE INTER THE DATE ?) ED ATE 

INS NTEXT '•'•WING SECT10N-TMJ a NFF"«.CfiLE>"l"«DATE>"l' 

ECHO APR SYM OFF 

ft EL MOD MOD 

SEL LAY 6 

NUM = 52 

CHECK = INC: * 100 
INS POi: < # 

REPEAT , _ 

X A FOX < NUM ) > Y 0 Z! F 0 X ( NUM - l ) 3 » ' # 

NI.IM = NUM + 100 
UNTIL (NUM .GT. CHECK) 

•C CR > 

ZOO DRA win: <VAR> 

ROT ENT COPY AX 2 70 : < VAR > 

ZOO DRA ALL 
ECH LAY 6 

ERASE ENT: PWIN < VAR > 

ECH LAY ALL 

^.C;||p|p s “M 11 ^ ^ 

READ (DO YOU WANT TO SURFACE THE WING ?) &SURF 
IF (& SURF . f: 0 . "N") END 

<# 

.( g 

<# ROUTINE TO SURFACE WING BEGINS HERE 

<# 

CONVERT ENT :BSPL LAY 0 VWIN NAME RIGHT 

SEL LAY 10 

NUM = 50 

POLE = 1 

CHECK = INC* 1 00 

TOTAL = 49 

#SURF 

REP GEN SFTiLE : XfFOX (POLE) 3YGFOY (POLE) 37-lFOX (NUM) >» <# 
XGFCiX (POLE * 100) > YI FOY (POLE+100) >ZG FOX (NUM+100) > 


f p ;> 

POLE = POLE + 1 
UNTIL (POLE .EQ. TOTAL) 
POLE = POLE + 52 
NUM = NUM + 100 


VAR 


TOTAL = TOTAL +100 

IF (NUM+300 . LT. CHECK) GOTO SURF 



WING GENERATOR - SAMPLE RUN 


<# 

<# 

•(# 

<# 

# 

FMHJ.H.MQP1 IS 


THE HARDFILE 


#03# RUN NEW P r NPF .WIND 


PLEASE ENTER THE WING SPAN IN INCHES 
WING SPAN *? 1 200 

PLEASE ENTER THE LENGTH OF THE ROOT CHORD IN INCHES 
root CHORD LENGTH ? 240 

PLEASE enter the tip chord length in INCHES 
tip chord length ? 120 


PLEASE ENTER THE SWEEP ANGLE IN DEGREES 
SWEEP ANGLE -? -30 


WHAT IS THE AIRFOIL SECTION FILENAME (DEFAULT = (P.THJ.FUlL) 
P.THJ.FSWEPT 

THERE ARE! 0 AIRFOILS CONTAINED IN THE EJ. LE. 

IlRFOlL SECTION (1) HAS BEEN UPLOADED. 

HOW MANY SPAN STATIONS WILL (1) APPEAR 

{CANS) 1 

PLEASE INDICATE THE # (1) SPAN POSITION FOR AIRFOIL 1 
SPAN POSITION =? 0 


F'LEASE ENTER THE TWIST ANGLE 


IN DEGREES FOR SAME 


TWIST ANGLE =? -5 


airfoil section 

AT HOW MANY SPAN 


?■) HAS BEEN UPLOADED. 
STATIONS WILL (2) APPEAR? 


{ CANS > 1 

PLEASE INDICATE THE * 111 SPAN POSITION FOR AIRFOIL Z 


SPAN POSITION -? 150 



PLEASE ENTER THE TWIST ANGLE IN DEGREES FOR SAME 
TWIST ANGLE =? -2.5 

AIRFOIL SECTION (3) HAS BEEN UPLOADED. 

AT HOW MANY SPAN STATIONS WILL (3) APPEAR? 

-C&ANSM 

PLEASE INDICATE THE # (1) SPAN POSITION FOR AIRFOIL 3 
SPAN POSITION =? 300 

PLEASE ENTER THE TWIST ANGLE IN DEGREES FOR SAME 
TWIST ANGLE = ? 0 

AIRFOIL SECTION (4) HAS BEEN UPLOADED. 

AT HOW MANY SPAN STATIONS WILL (4) APPEAR? 

T&ANSM 

PLEASE INDICATE THE # (1) SPAN POSITION FOR AIRFOIL 4 
SPAN POSITION =? 450 

PLEASE ENTER THE TWIST ANGLE IN DEGREES FOR SAME 
TWIST ANGLE -? 2.5 

AIRFOIL SECTION (5) HAS BEEN UPLOADED. 

AT HOW MANY SPAN STATIONS WILL (5) APPEAR? 
i &ANS > 1 

PLEASE INDICATE THE. # <1) SPAN POSITION FOR AIRFOIL 5 
SPAN POSITION =? S00 

PLEASE ENTER THE TWIST ANGLE IN DEGREES FOR SAME 
TWIST ANGLE =? 5 

DO YOU WANT TO SEE THE AIRFOIL COORDINATES (N) 

DO YOU WANT TO SEE THE POINTS PLOTTED (N) 

PLEASE ENTER THE DATE () 

1/10/90 

DO YOU WANT TO SURFACE THE WING IN) 

# 03# < V A R > D Y N a ir. x c VI Ew SPEED 5: view d 
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#03#<VAR> 

#03#<VAR>DYNamic VIEW SPEED 5: 
#03#< VAR > 

#03# 

#03#SEL LAY 1 

#03# INS LAB "UPPER BOUNDARY" ’• 
#03#< VAR > 

#03#SEL. LAY 5 

#03# I NS BSF‘1 TAG ^ LOW: MODEL 
#03# I NS LAB "LOWER BOUNDARY": 
#03# 2 00 in DRA ALL 


v i s w d 


DRAW./ MODEL ent NEar 


loc: NEar <VAI\>ddddd 
DRAW/MODEL ent NEar 


< VAR>d DRAW loc 


< VAR>d DRAW loc 


#03#ERAse ENT : MODEL 
#03#ERAse ENT : MODEL 
#03#ERAse ENT : MODEL 
#03#ECHo TAG OFF 


ent Vie NAME TOP TA.<? HI LOW 

ent Vie NAME ISO TAs HI LOW GUARCH 

ent Vie NAME RIGHT TAs GUARCH 


#03#SEL MOD Dr a 

S' e 1 p c t e d mode i s D R A W . 
#03#ECHo APR SYm OFf 
#03#SEL MOD Mod 


Selected mode is MODEL. 
#03#SEL LAY 6 


#03#ROT 


#03# . , 

#03# ZOO DRA WIN: DRAW loc < VAR .■>- \ 

PM 'I COPY A X 2. 7 0 : MODEL ent <VAR>d 
MODEL 1 oc 
COPY AX 770: 

MOD El. ent 
MODEL ent 
MODEL ent 


MODEL loc 

0 

POi 

di 

# # # # # 0 3# •( VAR 

>ROT ENT 

MODEL loc 

0 

POi 

d : 

MODEL loc 

0 

RO i 

d : 

MODEL loc 

0 

POi 

d : 

MODEL loc 

o 

POi 

d 

#03# < V A R > 




#03# ZOO DRA 

AL 

L 



DRAW loc DRAW loc 


> #03#< VAR> = ####ROT ENT 
MODEL ent d 
d 
d 
d 


dd DRAW 1 
COPY AX 270 


#03# EC H LAY 6 

#03#ERAse ENT*. MODEL ent PWin DRAW loc <VAR ,»dddddd 
#03#ECH LAY ALL 


#03# DO HARd 
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The computer code asks the user to input severa 


pieces of data before the weight iteration begins. 


1 . Lift coefficient: 


The value of CL asked for here 
is the level flight at 
design altitude. 


2. Mach Number Range: 


This input is the values of the 
lower and upper mach number 
range as well as the step 
interval. 


3 . payload weight: 


The weight of payload to be 
carried in pounds. 


4. Thickness ratio: 

5. Sweep angle: 

6. Limit load factor: 

7. Aspect ratio: 

8. No. of propellers: 

9. Number of blades: 


The maximum thickness to chord 
ratio entered as a decimal 
value . 

The wing quarter chord sweep in 
radians . 

The value where yielding will 
occur . 

Value for main wing. Wing 
weight estimates may be low for 
high aspect ratios. 

Total number of propellers. 

Number of blades for each 
propeller. 


10. Propeller diameter: Maximum propeller diameter in 

feet. 


After this data is entered, the program sets values 
for constants at the design altitude. These initial 
values are given in Figure A. 3.1. The program will then 
iterate the component weights over the range of Mach 
numbers specified. A flowchart of the iteration method 


is given in Figure A. 3. 2. 

The iteration method is based on a graph of 
calculated gross weight (WG) versus an initial weight 
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. * o o Tti 0 4 5 line 

guess (WI) as shown in Figure • . 

• * i s np where the initial 

starting at the origin is the line whe 

q uess equals the iterated weight. At low initial 

weights, the calculated weight will be greater. As 

initial guesses increase, the calculated weight 

decreases and the line traced out will cross the line 

where the initial and calculated weights are equal 

This is the point where the weights have converge . 

4 - tnii-ial weights (DW) until this 
program will increment initial g 

At this point the initial weight will 

crossover occurs . 

be r eset to the last value (») where the calculated 
weight is greater than the initial weig 
increment will be reduced and the procedure will begin 
again until the calculated value is within live pounds 

of the initial value. 

The iterations begin with an initial weight guess 

. ■ u is the payload weight. The lift 
of 1000 pounds, which is p 

to drag ratio tor the initial sizing is calculated 
lir st. This calculation is accomplished as tollows. 
total drag is the sum ot pressure and sKin friction 
drags on the wing, tail and fuselage plus the induced 

drag on the wing. 


= C + K'C * + K"C l 


(A. 3.1) 


where : 


K' 
K' ' 


i / ( ir * e * AR ) 

Average value of NACA 2412 
obtained from Reference 1 
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The skin friction drags, C , on the wing, tail and 
fuselage are determined based on flat plate, turbulent 
flow, and are referenced to the wing area by the 
following formula: 


c = [.074 / n b ; j ] s wet /s ref 


(A. 3. 2) 


The total skin friction drag is then multiplied by 
1.25 to account for pressure drag and mutual 
interference as suggested on p. 2-14 of Reference 24. 

The sum of the skin friction drags is the drag 
coefficient at zero lift, C.. The induced drag is 
determined using a span efficiency factor, e, of .85, 
and viscous drag is calculated based on an estimate from 

_ _ i ail three are summed to obtain 

p. 479 of Reference 1. All rnree ««= 

the total drag coefficient, C*. 

The power required is calculated based on the 

current weight estimate. This value is simply the 
velocity times the drag force. The power required 
is multiplied by a factor of 1.22 to account for 
equipment efficiency, motor (97%), DC conversion 
(98%), propeller (85%) and 10% extra power used for 

climbing . 


PR = 1 . 22 (WI/LD) ( V) 


(A. 3. 3) 


All of the component weights are now evaluated 
and the size of the external rectenna, if required, is 
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determined. The rectenna sizing is determined by the 
following methodology. The total rectenna area required 
to absorb the necessary power is calculated as: 

AN = PR / .0872 (A. 3.4) 


where : 


AN - Total Rectenna Area (ft 2 ) 

PR - Power Required (horsepower) 
.0872 - Power Density (hp/ft 2 ) from 
Reference 6 


If this area is less than or equal to the wing area 
then no external rectenna is needed. The weight of the 
rectenna on the bottom of the wing surface is estimated 
as .076 pounds per square foot also suggested by W. 

Brown in Reference 6. This weight includes the weight 
of the thin film covering and the reflecting plane o 
the rectenna. If the required rectenna area is larger 
than the wing area, the difference is the external 
rectenna area. The external rectenna is assumed to be a 
flat disk. The diameter is determined and the power 
required is re-calculated to include the drag of the 


external rectenna. 

When convergence is reached, the program will 
output data to the monitor, printer and the specified 
data file. The output to the monitor and printer are 
broken down into component weights, span, lift-to-drag 
ratio, power available, wing area, wing loading and 
ultimate load factor. The output to the data file is in 
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ASCII format. The 
so that each time 


data file is opened as an append file 
the same filename is used the new data 


will be appended to the original file 


This output 


includes the Mach number, lift coefficient, gross 
weight, wing area, external rectenna area, lift to drag 
ratio, optimum lift-to-drag ratio, span and power 


available . 
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Weight Estimation Program for HALE Aircraft 

Written By: Scott B. Sandler 

Date: Revised November 18, 1989 V2.1 


List of Variables and Constants 

CL - Coefficient of Lift 

M - Mach Number 

LD - Lift to Drag Ratio 

AR - Aspect Ratio 

TC - Thickness Ratio 

SW - Sweep Angle (rad) 

WP - Payload Weight (lbs) 

CD - Drag Coefficient 

WA - Wgt/Thin Film on Wing (lbs) 

_ _ . _ ^ a / *1 


NRE - Reynolds Number 
S - Wing Area (ft2) 

MU - Viscosity 
B - Span (ft) 

:BAR - Mean Cord (ft) 

N - Load Factor 
CO - Drag at Zero Lift 
NP - Number of Propellers 
AN - Total Rectenna Area (ft2) 


WA - wgt/inin rum uu /ihc^ 

AX - Area External Rectenna (ft2) WX - Wgt E £e™al Rectenna^lbs^ 
DR - Diameter/External Rectenna CR - Skin Friction p u „ laae 
CW - Skin Friction on Wing CU - Skin Fri cti°n g 

F$ - Output Filename OFT “ Optimum L/ 

* ^ . — ■ , , cu _ UAri Tail Area 


190 REM SV - Vertical Tail Area 

200 REM BV - Vertical Tail Span 

210 REM CH - Skin Friction on HT 

220 REM F - Iteration Flag 

230 REM 

240 REM *** define constants *** 


SH - Horizontal Tail Area 
BH - Horizontal Tail Span 
CV - Skin Friction on VT 


270 RO = 3.211E-05 :REM denisty at altitude (lb/ft3) 

non i/tt sz ^ 1 5E-07 : REM Viscosity 

290 SS = 1003 2 :REM local sonic speed (ft/s) 

300 I - 85 :REM Efficiency Factor 

iio cm -35 : REM Gust velocity (ft/s) 

320 AU = 500 :REM Uninstalled Avionics Weight (lbs) 

_ q :REM taper ratio „ . . 

”2 m- 007 :REM Viscous Drag Coefficient estimated from a 

340 KV - . 00 / naca 2214 airfoil 

„ n »p H _ = 2 :REM Horizontal Tail Aspect Ratio 

360 ARV = l i 33 :REM Vertical Tail Aspect Ratio 

365 ARC = 4.7 „ ^ , 

970 F - 1 :REM Set Fla< ? to 1 , , . 

380 L = 70 : D = 4 :REM Length and Diameter of fuselage (ft) 

390 INPUT "Output Data File Name [drive ]:[ path\f ilename ]"; F$ 

400 OPEN F$ FOR APPEND AS 1 LEN=2000 
410 REM 

420 REM *** data entry *** 

430 REM 

450 INPUT "Lowery Upper Mach Number and Increment (ML , MU , I ) " ;M1 ,M2 , I 

460 INPUT ’’Payload Weight (lbs) 11 ;WP 

470 INPUT "Thickness Ratio" ?TC 

480 INPUT "Sweep Angle (radians )" ?SW 

490 INPUT "Limit Load Factor" ;N 

500 INPUT "Aspect Ratio" ;AR 

520 INPUT "Number of Propellers" ;NP 

530 INPUT "Number of Blades per Propeller ;NB 

540 INPUT "Diameter of Propeller" ; DP . 

550 DW = WP : WG=WP : WI = WP : REM weight increment 

555 CLS 
560 REM 

570 REM *** begin iteration *** 

580 REM 


585 PRINT "Iterating Mach Number: 

590 TO R M = Ml TO M2 STEP X 

|oo : REM wing loading 

610 GOSUB 1100 : REM Calculate t div erged 

Ifo ll lisU-T)-< 5 THEN 750 ': REM check for convergence 

H° 0 - DELTA 

III IF DELTA = ODELTA THEN LW=WI ELSE 710 

670 F = 2 

680 WI = WI + DW 

690 ODELTA = DELTA 

700 GOTO 600 

710 DW = DW/2 : F - 1 

720 WI = LW 

730 DELTA = ODELTA 

740 GOTO 600 

750 GOSUB 810 wr-VT M 

760 WI = WP : DW = WP : F-l • NEXT^M 

770 PRINT "Iteration complete 

780 CLOSE 1 
790 STOP 

800 REM . 

810 REM *** print routine *** 

III S-.Tn.M,C l ,WO:S ! AX,LD;OPT = -T/C> 

840 LPRINT "Component Weights ar 

850 LPRINT " ; WW 

" ;WA 
" ?WX 

11 ;WT 

" ; WV 
" ;WC 
" ; WF 
" ;WL 
" ;WH 
";WE 


O DU . 

860 LPRINT "Wing Weight: 

870 LPRINT "Weight Thin Film. 

R80 LPRINT "Weight Ext Rectenna. 

IIS LPRINT "Horizontal Tail Weight: 

ISS LPRINT "Vertical Tail Weight: 

910 LPRINT "Controls Weight: 

920 LPRINT "Fuselage "W*: 

930 LPRINT "Landing Gear ” ei( 3^: 

940 LPRINT "Hydraulic System Weight. 

Ill LPRINT "Electrical System Weight ^ 

960 LPRINT "Motor Weight (Inst & In f t *, ) . WT ^ WM 
970 LPRINT" Gearbox Weight • 1 

980 LPRINT "Propeller Weight: 

990 LPRINT "Payload Weight: 

1000 LPRINT "Gross Weight: 

1010 LPRINT "Initial Guess: 

1020 LPRINT "Wing Area: ?S, 9-^; ft." 

1030 LPRINT "External Rectenna Area. . 

1040 LPRINT "Wing Loading: ,WS 

1050 LPRINT "L/D optimum ;OPT 
1060 LPRINT "L/D Ratio: ;LD 

1070 LPRINT "Power Reguired. ,PR 
1080 LPRINT: LPRINT 
1090 RETURN 

1110 REM *** weight calculations *** 

1120 REM 

\lll GOSU B ° 1470 : REM “Wi?5) ift t0 dra9 rati ° 
lll° 0 g I ( ({k*!- N 65) * " I ( AR/COS ( SW ) ) * . 57 ) * (S/100) 


: " ; WTR 
"?WB 
" ; WP 
" ; WG 
";WI 

ft. 


. 61 ) 



1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 


ww 

WT 

WV 

WC 

WF 

WL 

WH 

WO 

WE 

PG 

WB 

AN 

AX 

IF 

WX 

WA 


= .75 * (96.948 * ( W jj * \ j ( TS*! . 3976* < ( ARH/ (12*TC )) * . 5 ) ) * . 458 ) 

: :]l : ! 98^5* |^K1' ^87 ) * ( ( • 01*SV) * 1 . 2 ) * ( ( ARV/ ( 12*TC) ) A . 5 ) ) ) 

: ^75**[ ( 200* ? ( ( (Kl) “ . 286 ) * ( ( L/10 ) A .857)*((2*D)/10)*1.36302) A 1.1) 

I ! 0005M WI *1.28) : REM hydraulic system weight 

= 2.117*(AU a .933) 

= 4 26* ( (WO/ 1000 ) A . 51 ) -niHai Power Guess 

: • RE ” prop " ei5ht 

; S’SiX *" a »ori=o n t,l tail 

.2044 ; 3 REM .2044 lb/ft2 (film + structure) 


= (WI/WS) * .076 : REM 


weight of thin film 
on entire wing 


and reflecting plane 


DR = 
GOTO 
WA = 


(4*AX/3. 14159) * . 
1380 

AN * .076 : REM 


weight of thin film and reflecting plane 
portion of wing needed 


on 


REM Span of Main Wing 


llli m : ^2 . 575* ( PG/NP ) A . 922 ) *NP 
1390 WTR = - 2 5 * PG 

1400 GOSUB 1470 __ M . ecn 

1410 PR -1-22* (WVLD) * H |J 12 70 7 j REM check for power 
lil°0 WG t Vl m i WT + WF +WE + WH +WP +WX + WA +WC+WL 
1450 RETURN 

i5?0 REM *** calc lift to drag ratio *** 

1480 REM 

1490 S = WI / WS 
1500 SV = .1*S : SH * *2 S 
1590 B « SQR( AR * S ) 

1600 CBAR = S / B 
1610 NRE « ( RO*SS*M*CBAR) / MU 
1620 CF®. 074/( (NRE) A -2) 

1630 CW = CF * 2 
1640 BH - SQR ( ARH * SH ) 

1650 BV = SQR ( ARV * SV ) 

1660 NRE = (RO*SS*M*(SH/BH)) / MU 
1670 CF * *074 / (NRE A .2) 

1680 CH = CF * ( 2*SH ) / ® 

1690 NRE = ( RO*SS*M*L) / MU 
1700 CF=.074/( (NRE)*. 2) 
i 7 in ru = CF * (L*3 . 14159*D ) /S 
i ? 20 NRE = ( RO*SS*M* ( SV/BV) ) / MU 
1730 CF = .074 / (NRE A .2) g 


required convergence 

+WB + WTR +WM 


: REM Span of 
: REM Span of 


Horizontal Tail 
Vertical Tail 


1740 CV = CF * 

1750 CR = 0 

1760 IF DR = 0 
1770 NRE - (RO 
1780 CF = .074 / ((NRE*. 2)) 

1790 CR = CF * (2*AX /S) 

1800 CO * 1.25 

1810 CD = CO + 

1840 RETURN 
1850 END 


THEN 1800 

* SS * M * DR ) / MU 


* .. pw + CR t CH *■ CV ) 

UC^2) / (3.14159 * AR * E ) ) + KV*CL A 2^ 


n 4. KV*CL a 2) 



Figure A. 3.1 

First Iteration Weight Code Constants 


Altitude 

Air density 

Viscosity 

Speed of sound 

Span efficiency factor 

Taper ratio 

Viscous drag coefficient 
Fuselage length 
Fuselage diameter 


100,000 ft. 3 

3.211(10 ) lb/ft 

3 . 15(10 ) slug/ ( ft s) 

1003.2 ft/s 
.85 
.9 

.007 
28 ft 
4 ft 




Convergence Method Flow Chart 

Figure A. 3. 2 
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The wing-body drag polar is determined from the 

. oa a brief overview is given 

methodology in Reference 24. A brie 

here. 


c = C„r. + K'C* + K"(C l - 
U D DO 


) 


(A. 4.1) 


where : 


'■'DO 

K' ' 

K' 

C LMIN 


- zero lift drag j e "v S f °a ci?ve 

body (determined from C L vs. 

_ viscous drag factor* 

- inviscid drag factor . # , 

- lift coefficient at minimum drag 

may he approximated by l*** attack 
coefficient at zero angle of attack 


The zero lift drag coefficient term is comprised of the 
zero lift drag on the wing and of the fuselage. They 
are evaluated as follows: 

/ t4,„ o , s (A. 4. 2) 

_ _ c*[l+L(t/c)+100(t/c) ]R s we t' S ref 

DO (wing) f 


where : 


c - turbulent flat plate ^ in friction 
£« - maximum thickness location factor 

t/c - maximum wing thickness factor 

r - lifting surface correlation factor 

c - wetted wing area 

fwet _ reference area (wing planform) 

^ref 


Similarly, the drag on the body is: 

c = C^[l+60/(l R /d) 3 +.0025(l B /d)) S s /S B 

^DOtbody) f D 

where: ^ _ length of the body 

d B - diameter of body 

c - wetted area of body 

! s - maximum cross sectional area 
B 


(A. 4. 3) 
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The value for the inviscid drag factor was 


also determined from Reference 24. 


where: 


1 

K' 

e*AR 


(A. 4. 4) 


e = e' [l-(d/b) * ] (A. 4. 5) 

d/b - body diameter to span ratio 

e' - Weissinger wing planform efficiency factor 

The value of K" was determined from Reference 30. 


After evaluation, the final equation for the wing- 


body drag polar reduces to : 

C D =.0154+. 02167 C^ + . 1 ( 2 (A. 4. 6) 
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Preliminary Design 


To begin the analysis, simple wing models 
consisting mainly of the wing box were used. The wing 
had a chord (c) of ten feet, a span (b) of 250 feet and 
an aspect ratio (AR) of twenty-five. To minimize the 
size of the computational model, only one - half of the 
wing was modelled due to symmetry. Translations in the 
X, Y, and Z directions were restrained at the aircraft 
centerline. Aerodynamic forces were represented by an 
elliptic load distribution of 5,000 pounds on each wing 
panel. The wing was constructed out of graphite - epoxy 
with a modulus of elasticity of 20 X 10* psi, a 
Poisson's ratio of 0.28 and a density of 2.91 slugs/ft 3 . 

After running this ANSYS model a maximum deflection 
of 2.9 feet, an average maximum shear stress of 0.271 ksi 
and an average equivalent stress of 114.6 ksi were 
obtained. The ultimate stress for graphite - epoxy is 
69.9 ksi. The vertical reaction forces, 4,690 pounds, 
differ from the wing loading, 5,000 pounds, by only six 
percent. This difference in loading is due to the fact 
that the pressure was placed perpendicularly to the 
surface of the upper panel, which is inclined relative 
to the wing refrence panel, hence the vertical force 
component differs by the cosine of the angle. 

These results indicate that the wing will fail due 
to the high equivalent stresses. Furthermore, the weight 


Of each wing panel is 5,308.7 pounds, which exceeds the 

wing weight allowance of 2700 lbs. 

Next, it was necessary to expand the structural 

model to include both leading and trailing edges, 
thereby making the model more realistic. Second, the 
thickness of the wing material was decreased in order to 
reduce to overall weight of the wing. The thickness of 
the first four sections out from the centerline was 0.1 
inches, of the middle four sections 0.01 inches, and of 
the last four sections 0.0025 inches. The same material 

properties were used in each section. 

The large decrease in material thickness, with a 
commensurately large decrease in cross sectional area, 
resulted in a dramatic increase in stresses every fourth 
section. For this configuration the stresses obtained 
were 11.8 ksi for the average maximum shear stress and 
206.25 ksi for the average maximum equivalent stress, as 

shown in Figures A. 5.1 - A. 5. 5. 

The stresses proved to be three times as much as 

the ultimate stress, 69.9 ksi. Also, the minimum 
thickness of graphite - epoxy is approximately 0.03 
inches as stated in Reference 28. 

To alleviate the above conditions the panel 
thickness was reduced gradually, which yielded a 
gradual change in the stress distribution, second, to 
reduce the maximum stresses flanges were added to the 
wing spars, and not just webs as previously shown. 
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Finally, in order to have a more realistic model, a more 
accurate elliptical loading distribution which correctly 
represents the equivalent bending moment, was input. 

This will result in the correct I beam thicknesses and 

cross sectional areas • 

The equation of the spanwise elliptical loading is 

w(y)= Wo(l-(y/(b/2))2))’, 


where Wo is the maximum load (found at the centerline), 
b is the span of the wing, y is the distance along the 
span and w(y) is the load in lbs/ft found at varying 
points along the wing. 

The bending moment is found by integrating the 
above equation twice. To accomplish this task a program 
was written to perform the integration, as shown at the 
end of this Appendix. To test the accuracy of the 
program the equation of elliptical loading was changed 
to an equation of uniform load and then later to an 
equation of load concentrated at one point. All three 
runs were found to be accurate, as shown in Figures 

A. 5. 6 - A. 5. 14. 
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Model /I 


After completion of this preliminary model and 
analysis, new wing geometry and weight was provided by 
the aerodynamics group. The aspect ratio (AR) was 15, 
the thickness to chord ratio was 0.12 and the total wing 

weight was 2704.644 pounds. 

Given the total wing area to be 3633.8 square feet 

and the minimum thickness of graphite - epoxy to be 
0.035 inches, the skin weight would be 2270 pounds. 

The skin weight left an unreasonably low weight of 434 
pounds available for use in interior structural support. 

Due to this high weight, it became necessary to 
turn to another material to construct the skin. Aluminum 
was chosen because it is possible to reduce the material 
thickness to 0.01 inches. But with aluminum the density 
rises to 5.42 slugs/ft’ yielding a skin weight of 1059 
pounds. This left 1645.6 pounds to be used for the 
interior structural support. Next, it was necessary to 
determine the respective sizes of the two I beams in the 
wing. The forward wing spar was sized using 67 percent 
of the load, and the rear spar was size using 33 
percent, as shown in Figures A. 5. 15 - A. 5. 20. For this 
initial analysis it was determined that the plane would 
be designed for an ultimate load factor of three, which 
later proved to be much lower than required. 

Given the weight of the plane, from the 
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aerodynamics group, to be 6718 pounds and a load factor 
of n =3, the load the wings were designed to 

ult 

withstand was 20,154 pounds. 

To get the maximum Wo one can use the relationship 

L = (ir*Wo*(b/2) )/2. 


Solving for Wo yields, 

Wo = (L*4)/(»*b) 

where L is the load of 20,154 lbs and b is the span of 
233 ft. Hence, the final result is a Wo of 110.13 lb/ft. 
As stated previously the front wing spar carries 67 
percent of the load and the rear wing spar carries 33 
percent of the load. Hence, one gets a Wo |fM|tj of 6.15 
lb/in and a Wo of 3.03 lb/in. 

9 (rear) 

By using the same integration program to integrate 
the elliptic load function, 

w(y) = Wo(l-(y/(b/2))2)^, 

one gets the moments on the front and rear wing spars. 
Additionally, from a Lotus 1-2-3 output one gets the 
average moment every 11.5 feet (in other words every 
span section) . 

To find the thickness of the front wing spar one 
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can use the following equation. 


s min = < M W/ (a) all 

in the above equation S min is the minimum section 
modulus, (M) max is the maximum moment and (o) all 1S the 
allowable tensile stress. Knowing that S=I/c, where c is 
the maximum distance from the neutral axis, and assuming 
h = b = 0 . 12 * (chord) one can obtain the equation for the 

moment of inertia , 

1^ = (b*h 3 /12) ” ( 2 * b i*^i / 12 ) 

The same procedure is used to find the thickness of the 
rear I beam. The results are given in Figures A. 5. 21 s 

A. 5. 22. 

Nonlinear deflections 

Using Euler - Bernoulli analysis from ANSYS one 
gets a vertical deflection of 198 inches. Now using 
Figure 2 on page 1053 of Reference 17 a deflection to 
length ratio d/1 can be calculated in order to compare 
the results to a nonlinear deflection. Hence, one gets 
d/1 equal to 0.14, K equal to 1.44 and t equal to 0.1 
inches. Here K is equivalent to 
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From the graph on page 1053 the Euler - Bernoulli 
theory shows a close relation to the Rhode's method in 
this case. However, one should note that equation 
(A. 5.1) j_s valid for a uniform load distribution and not 
for an elliptical distribution, hence the validity of 
its use in this case is uncertain. Also, one should 
note that the Rhode's method of solving nonlinear 
deflections was the only one practical. 
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Preliminary Loads Analysis: 

The first step in finding the correct load 
distribution along the wing is to find the proper 
airfoil section lift coefficients along the span of the 
wing - Given that the lift over the wing is elliptical, 

one can use the equation W(y) - Wo(l-y/<b/2) )) - 

AlEO given that Wfy) = CL* 0 . 5 *den.v.v.c (y) *y, one can 

equate the above equations to solve for Cl. The 
following conditions were used, a Wo of 110.14 lb/ft, 
semi-span of 115 ft, a density of 3.211B-5 slugs/ft . a 
velocity of 603.2 ft/s, a mean chord c(y) of 15 ft and a 

of 11 5 feet. The final equation will then 
span section of 11.5 ieei-. 

be 




2*(Wo(l~(y/(b/2) ) A 2 ) ) 


CL = 


(A. 5. 2) 


( den* V*V*c ( y ) ) 


As one moves out from the root to the tip of the 

• • +• with Wo, and Wo varies 

wing the lift coefficient varies witn « 

elliptically. 

A fter finding Cl along 10 evenly spaced positions 
along the span the lift coefficients were given to the 
aerodynamics group. The aerodynamics group then 
calculated the corresponding pressure coefficients Cp 
chordwise along each of the cl's given. 
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However, this analysis (load = 3W) led to a Cl of 
1.8 (the aerodynamic group's CL = 0.6 for L-W). Knowing 
that the plane would stall at about a CL of 1.0 the 
situation had to be reevaluated. It was decided to solve 
for chordwise pressure coefficients across the span. 
Hence, one had to assume that the maximum CL encountered 
during gust conditions would be a CL = 1.0. But this 
would give a different velocity than that used by the 
aerodynamics group (M=0.44). Hence, to get the correct 
wing loading of 20154 lbs a new velocity was required, 
which turned out to be 603.2 ft/s. Knowing the velocity 
one can then calculate q (dynamic pressure) and Cl along 
the span to determine the pressures on the wing. The 
corresponding results are given in Figure A. 5. 24. 

Using the average panel pressures the results were 
entered into the ANSYS model with a gravity load, the 
results were as follows: vertical deflection of 14 feet, 
maximum equivalent stress of 30 ksi and maximum shear of 
733 psi. One analysis which has yet to be considered is 
the effects due to torsion. Since ANSYS has no direct 
commands to check for torsion it was determined if the 
torsion in wing was too high by performing a twist 
analysis . 

Another note of interest is that the deflection in 
the chordwise direction was 7.33 inches, which is 
believed to be too high for such a model. This 
deflection could be due to the low stiffening in the 
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model . one way to compensate for the lac* o£ stiMenmc, 
in the chordwise direction was to replace the bar 
elements along the leading and trailing edges with beam 
element, since beam elements can better represent the 
stiftening components along the desired directions, 
beams along the leading and trailing edges were made 
higger (using the average real constants from the front 
end rear I beams respectively, , which increased the 
weight to 2319 lbs. 

The bar and beam elements along the leading and 
trailing edges were checKed for twist. The desired 
angle of twist was plotted versus the semi-span of the 
wing, m order to determine the angle one must tabe the 
arcsine of the change in chordwise direction divided it 
by the chord. Hence one gets 

Sin(x) = delY/chord (A. 5. 3) 


where delY is the change in displacement from the 

trailing to leading edge. 

one can readily see from Figures A.5.25 S A.5.26 

that the angle of twist reaches a maximum of 0.17 

4 -Vial- the winq is experiencing very 
degrees. This means that the wing 

i. 4 -V 10 static torsion on the wing is 

little twist and hence the stati 

also very small. 
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Second Model 


Since the first model seemed to lack sufficient 
stiffness a second model was developed. The geometry 
will be maintained from the previous model but with an 
additional I beam placed in the middle of the wing. 
Again, the I beams will vary in thickness as shown in 
Figures 6.1.13-6.1.16. The old I beams (wing spars) 
had unrealistic dimensions for their flanges, which were 
21.6 and 14.4 inches wide. The flanges were redesigned 
and reduced to a width of 4.32 inches in the front, 3.74 
inches at the middle, and 3.17 inches in the rear, but 
the cross sectional area was kept constant, as shown in 
figure A. 5. 27. There will also be twice as many span 
stations, 5.75 feet apart. Finally, updated pressure 
distributions supplied by the aerodynamics group have 
been used to calculate the applied load distribution on 

the wing. 
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Loads Analysis II 


From the aerodynamics group's analysis of the 
chordwise pressure distribution one can establish a 
better chord and spanwise pressure distribution per 
panel on the new wing (160 surface panels). 

The chordwise pressure coefficient distribution for 
a CL of 0.743913 was provided by the aerodynamics group. 
The wing structural model had four chordwise panels, so 
the pressure distribution was approximated as shown in 

Figure 6.1.2. 

To determine if these are the correct pressures one 
can compare the calculated to the expected force and see 
how they relate. The comparison between the expected 
and calculated results are shown in Figure A. 5. 28. 
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Second Model 


As mentioned in Reference 28, graphite epoxy will 
not bond directly to aluminum. Graphite - epoxy will 
need an E glass filling between all joints where it is 
connected to aluminum. This will raise the weight of 
the wings significantly beyond the designated 2700 lbs. 

The wing weight of the new model with aluminum skin 
and no E glass filling is 2800 lbs. The weight is 
already 100 lbs over the designated allowance without 
the weight of the E glass taken into account. 

Further investigation on graphite - epoxy wing 
construction indicated that building a wing skin of 0.01 
inches is possible. Graphite - epoxy can be manufactured 
to about 0.02 inches, so, as stated in Reference 26, it 
might be possible to produce graphite - epoxy to 0.01 
inches in the near future. Hence, with this 
consideration in mind and the fact that this is still a 
first weight estimation, the skin was changed from 
aluminum to graphite - epoxy, reducing the wing weight 

to 2319 lbs. 
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Gravity Loading Analysis 


With gravity alone, one gets the following results: 
a maximum vertical deflection of -1.45 feet, an average 
maximum equivalent stress of 5183 psi and a maximum 
horizontal displacement of -0.81 inches, as shown in 
Figures A. 5. 29 - A. 5. 34. 

Pressure Loading Analysis 

With just a pressure loading, the average maximum 
shear stress is -2.0 ksi, the average maximum 
equivalent stress is 39.7 ksi, the maximum vertical 
displacement is 10.9 feet, and the maximum horizontal 
(x) displacement is 6 inches, as shown in Figures A. 5. 35 
- A. 5. 44. Comparing the results to those taken on the 
model which had only 10 span stations, larger flanges in 
the I beams and no middle I beam, the results were as 
expected. Since there was more interior stiffening in 
the model, the deflections were reduced from 16.5 to 
10.9 feet (vertical direction) and from 7.33 to 6 inches 
(horizontal direction). The shear stresses were reduced 
from 2184 to 2005 psi. However, the average maximum 
equivalent stress increased from 33,153 to 42,138 psi. 
This increase is attributed to the width reduction of 
the I beam flanges from 21.6 to 4.32 inches. But, since 
graphite - epoxy's ultimate principle stress is 69.9 
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ksi , the wing is still acceptable. 
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1 REM Program #1 

2 REM This program allows you to take the double 
10 REM Elliptical load distribution 

20 OPEN "a: amoment.dat" FOR OUTPUT AS 1 LEN-2000 

30 INPUT "what is the Max X value"; B 

10 INPUT "how many sections"; NTOT 

50 DIM LOD (NTOT- 1) , Y(NTOT-l) , SHR (NTOT- 1) , MOM {NTOT- 

SO DELY = B/NTOT 

70 LOD { I ) = 0 

30 Y(I) =0 

90 X=0 : GOSUB 1 000 : LOD ( 1 } =Y 
100 FOR 1=1 TO NTOT 
110 Y(l)= ( I - 1 ) ^ DELY 
-120 X - Y { I ) 


integral 


1 ) 


of 


an 


130 GOSUB 1000 

140 LOD ( I ) = Y 

150 NEXT I 

160 Y { NTOT-l ) =B 

il V 0 LOD { NTOT* 1 ) = 0 

180 REM compute shear 

‘190 SHR (NTOT-1) =0 

200 FOR K-NTOT TO 1 STEP -1 

.210 LAYE= (LOD(K)-LOD(K-l) ) /2 

£2C SHR (K)=SHR(K-l)-( LAVE * DELY ) 

.230 NEXT K 

*240 REM compute moment 

250 MOM (NTOT- 1 ) =0 

$26 0 FOR K - NTOT TO 1 STEP -1 

1270 MAYE= { SHR (K ) - SHR ( K - 1 ) ) / 2 

280 MOM (K) =MOM ( K - 1 ) - (MAVE* DELY ) 

290 NEXT k 

300 FOR X = NTOT- 1 TO 1 STEP -1 

310 PR I NT# 1 , Y ( X ) ; LOD ( X ) ; SHR { X } ; MOM ( X ) 

320 NEXT X 
330 END 

1000 IF X > B THEN X^B 

1001 Y = 6.15M <i-(X/B) '2) ' (1/2) ) 

orx, DV 1 the equ3Zlon you Wish to integrate in the above 


line 


equation 


1 REM Program #1 

2 REM This program allovs vou to take the double mtprtrai 

10 REM Elliptical load distribution integral of an equation 

20 OPEN' "a : amomen t . da t '' FOR OUTPUT AS 1 LEN = 2000 
30 INPUT "what is the Max X value"; B 
JO INPUT "how many sections"; NTOT 

50 DIM LOD (NTOT-1 ) , V (NTOT-1 ) , SHR (NTOT-1 ) , MOM (NTOT-1 ) 

50 DELY=B/NTOT 11 

70 LOD ( I ) = 0 
30 Y(I)=0 

90 X = 0 rGOSUB 1000: LOD ( 1 ) = Y 
100 FOR 1=1 TO NTOT 
110 Y (I ) = (I- 1 ) *DELY 
120 X = Y ( I ) 

130 GOSUB 1000 
110 LOD ( I ) =Y 
150 NEXT I 
160 Y (NTOT-1 > = B 
117 0 LOD (NTOT-1 ) =0 
180 REM compute shear 
‘190 SHR (NTOT-1 ) =0 
•200 FOR K- NTOT TO 1 STEP -1 
,210 LAVE* (LOD (K) -LOD ( K - 1 ) ) /2 
*220 SHR (K) = SHR ( K- 1 ) - ( LAVE' DELY) 

.2 30 NEXT K 

*2-10 REM compute moment 

.250 MOM (NTOT- 1 ) =0 

)260 FOR K=NTOT TO 1 STEP -1 

•270 MAVE* (SHR (K) -SHR (K-l) ) /2 

280 MOM(K) *MOM(K-l ) - (MAVE* DELY) 

290 NEXT K 

300 FOR X = NTOT- 1 TO 1 STEP -1 

310 PR I NT# 1 , Y (X) ; LOD ( X ) ; SHR (X) ;MOM(X) 

320 NEXT X 
330 END 

1000 IF X>B THEN X 

1001 Y=6. 15 


B 


( (1- (X/B) ‘ 2 ) ‘ (1/2) ) 
llll RETURN 1 the equaC1 ° n you Klsh to integrate in the above line 
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Figure A. 5. 2 
Shear Stress 



AN8YS 4.3A2 


Figure A • 5 • 3 
Shear Stress 
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Equivalent Principle Stress 
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Figure A. 5. 5 

Equivalent Principle Stress 
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Figure A. 5 . 6 

GRAPH#6.1.1 - (LOAD VS. LENGTH) 

uniform load 
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GRAPH#6.1.11 - (SHAER VS SEMI-SPAN) 

SSX IN REAR I BEAM 
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(MOMENT VS SEMI-SPAN) 
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GRAPH#6.1.13 - (LOAD VS. SEMI-SPAN) 

66* IN FRONT I BEAM 
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GRAPH#6.1.14 - (SHEAR VS. SEMI-SPAN) 

66* IN FRONT I BEAM . 
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Figure A. 5.21 

Front & Rear Spar Properties 


Section 

Front 

A ve .moment 
I beam 66% 
(lb in) 

Front 

Thickness 

(in) 

Rear 

Ave. moment 
I beam 33% 
(lb in) 

Rear 

Thickness 

(in) 

1 

3,595,672 

0 . 1 

1,771,526 

0 .11 

2 

2,768,463 

0 .08 

1,362,724 

0 .087 

3 

2,054,949 

0.061 

1,012,438 

0 . 07 

4 

1 ,460,337 

0.045 

719,483 

0 .047 

5 

978 , 358 

0.04 

482,020 

0.04 

6 

603 ,805 

0.04 

297 ,484 

0.04 

7 

329,770 

0 . 04 

162,472 

0.04 

8 

147 , 242 

0.04 

72,543 

0.04 

9 

44 , 366 

0.04 

21 ,858 

0.04 

10 

4,672 

0.04 

2,302 

0 . 04 


F ront 


Section 

1 

2 

3 

A 

** 

5 

6 

7 

8 

9 

10 


Figure A. 5. 22 
& Rear Spar Properties 


Front 

Front 

Area 

lx X 

( in 2 ) 

(in 4 ) 

2.16 

1 . 8E-3 

1.728 

9.216E-4 

1 . 318 

4 . 086E-4 

0.972 

1 . 64E-4 

0 .861 

1 . 152E-4 

0 . 864 

1.152E-4 

0 . 864 

1.152E-4 

0.864 

1.152E-4 

0.864 

1 . 152E-4 

0 . 864 

1.152E-4 


Rear 

Rear 

Area 

lx X 

( in 2 ) 

(in 4 ) 

1.684 

1.597E-3 

1.253 

7 . 90E-4 

1.01 

4.116E-4 

0.677 

1.25E-4 

0.576 

7 . 68E-5 

0.576 

7 . 68E-5 

0.576 

7 .68E-5 

0.576 

7 .68E-5 

0.576 

7.68E-5 

0.576 

7 .68E-5 



Pressure 


Figure A. 5. 23 

distribution over the wing 


Section 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


v; (y) 
(lb/in) 
9.167 
9.077 
8 . 894 
8.613 
8 . 223 
7.708 
7 . 040 
6.170 
4.993 
3.187 


Total 
Pressure 
(psi) 
0.0490 
0.0485 
0 .0476 
0.04605 
0.0440 
0 . 0412 
0.0376 
0.0330 
0.0267 
0 .017 


Top 

Pressure 
(psi ) 
0.3675 
0.0364 
0.0357 
0.0344 
0.0330 
0.0309 
0.0282 
0.0248 
0.0200 
0.0128 


Bottom 

Pressure 

(psi ) 
0.01225 
0.0121 
0.0119 
0.0115 
0.0110 
0.0103 
0.0094 
0 . 0083 
0.0068 
0.0042 



Figure 4.5.2- 
Pressure Calculations 
for Si Be 


5 paL lit: coel CL=0.6 L=2015^ 


, Six: 

t:c pent 


£.75 2.353117 

2 108.8835 

17 v ' 

i ' . 4> - ‘ % ’ 

3 106.6225 

28.75 2.281245 

4 103.1642 

40.25 2.207042 

5 98.34923 

51.75 2.104032 

6 51.37668 

63.2: 1.567701 

7 52.65:55 

]i.I5 1.750454 

8 72.8*41: 

86.25 1.558252 

c ic 

57.75 1.241133 

10 34.38575 

105.25 0.725620 


lift coef 

1.2716^6 

i . 2 L c 6 3 5 

1.232805 

1.192101 

1.137:40 

1.063365 

0.967:73 

0.842:65 

C .670721 

0.357565 


4 

8 C 

Cf 

Cp Cp 

2.2 

0.952 -0.524 

2.153 

C.946 -0.0227 

2.057 

0.932 -0.0243 

1.916 

0.512 -0.0252 

1.721 

0.683 -0.526 

1.478 

0.845 -0.0278 

1.188 

0.756 -0.0303 

0.541 

0.733 -0.0342 

0.44 

0.646 -0.0404 

-0.036 

0.505 -0.0526 


B 

E 

F 

Cp 

Cp 

Cp 

0.368 

0.27 

0. J»iE 

0.387 

0.244 

* ^ ' 

0.386 

0.255 

0.374 

0.383 

0.247 

0.357 

0.375 

0.121 

0.422 

0.374 

0.20S 

0.443 

0.36E 

0.181 

0.446 

0.361 

0.1*: 

0.415 

0.351 

0.052 

0.311 

0.336 

P.C1 : 6 

*0.015: 


Hi'.iii 1.65662. 


1 C. 

2 0 . 

■ 0 


: r. 

6 0 

1 0 

8 2. 

} V. 

10 "2 


id 

.4 1 
0632*6 
057343 
083-45 

,06582; 

.0:555" 

.048193 


:n 

::d t 

r •• c : 1 


c : 


Tt; fc?w i«fir Bottoi 

lux C Heir B Hi« l ? T5 - 
-C-.DJ557 0.315735 0.010553 0.014117 0.054.-. 5 
.831277 -0.00056 0.015653 0.005555 0.31448. 0-3;327. 
.C2;t05 -0.00055 C. 015655 0.01050: 0.31::-": 0.0;2-5-. 
r. f 'ipc .i- jjijj 0.015537 C. 010015 0.1 :6:-: 0. 0:0:5: 
,'mW1 -0*00*105 0.015374 0.005370 0.017115 0.047905 
*03*200 -0.00112 0.015171 5.038474 0.017571 0.0443,2 

-0,00122 0,014325 0.007342 0.018173 0.035834 

-0.0313? O.OI46-: 0.005501 0.016597 0.033861 

•'.526207 -0.00163 0.014238 0.003232 t.KKU 0.025735 

1 520649 -0.00224 0.013630 0.032308 -2.03073 0.912;*! 
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Figure A. 5. 28 

Farce Calculations 
for ii :£ 


•'rr 


i’v I' 


SECr.OI FECIT : nice ^ 

: 1A0.0S70 204.7426 

2 139.7361 204.2236 

3 129.3337 :c 2 .:::: 
if 127. 5659 2C-1 . i 

5 136.5377 153.5551 

6 124.7280 196.S1C1 

7 122. 5227 192.6885 

8 125.505? 185.5608 

9 126. 6“5: 185.3836 

10 122.3131 180.2265 

11 113.2636 17- . 2112 

12 11,. 6-85 167.5632 

13 105.3896 153.87/1 
” ; u 102.3310 151.1053 
15 5c .£1-396 141.0603 
* ( ? : : v> £ i J C 4 ? 1 1 


03 

:;i.i 6 s 6 

140.611? 


TO? 
i£AE 4 


43.1037. 
tj. 95574 
42.77920 
42 .45182 
42.01160 


BOTTOF 
Erf.? 5 

t ; 47 "jC 


jo J 

53.74468 

52.47375 

52.06477 


40.77652 

39.97070 

23.02925 

37.94251 

36.69713 


17 '3.13052 .1 5 . ■■‘•25 
15 67.84056 3?.;5i58 
15 52. 1--" 7 ;.8:3'5 


137.588? 

135.7642 
132.5431 
130.9040 
127.8212 
124.2617 
120.1831 
115.530* 
no. 2310 32.65834 
104.1863 31.81262 
37.25735 23.63651 
85.23647 27.24835 
?u::*5 24.26* 33 
66.36241 20.87401 
52.65469 16.38311 


i' ; .t5t?7 51.81846 
50.37066 
45.56337 
48.78673 
47.42814 
45.87141 
35.27647 44.09555 
42.07232 
35.76:78 
37.12114 
34 ,060tE 
30.4:866 
26.05252 
20.47885 


BOTTOM BCTTSI* 

mi 6 m: 7 

53.8756; 8.620712 
53.74166 8.595143 
53.47275 8.555631 
53.064:7 8.450364 
52.51451 8.402321 
51.81846 8.230554 
50.37066 8.155:^5 
45.96237 7.551140 
48.78673 7.815878 
47. 4281 1 - 7.588502 
45.87141 7.333126 
14.05355 7.055254 
42.01252 6.7: :668 
39.76576 6.362525 
37.12114 5.533382 
34.06046 5.145638 
30.45866 4.873266 
26.052:2 1.1 -503 
20.47889 3.276622 


BCTTOH CALCULATED 
FKOST 8 FOFCE 
-4.31037 642.167? 
-4.23557 641.561' 
-1.27750 635.32)7 
-1.24518 635.4706 
-4.20116 629.9226 
-4.14547 622.6257 
-4.07765 613.5506 
-3.95707 60z.56ij 
-3.90253 585.5621 
-3.79425 574.3946 
-3.66571 55:. 5656 
-2.52761 536.7375 
•3.36582 513.6678 
•3.16126 487.2127 
•2.96965 455. 'ill 
-2.72163 121.3138 
-2.13664 375.4581 
-l.CB 1 !"’ 328.5C.iO 
-1.63E21 262.6588 


m • C902 21.3772? 5.58084? 


1 r.uuc .r 


EI2ECTE2 

FCF.CE 

643. 661 8 

642.0436 

638.6131 

623.3::: 

627.2536 

615.0384 

608.5102 

556.B?’? 

582.8205 

566.5522' 

54? .5922 

526.7787 

502.615- 

475.0526 

443 .i555 

406.6365 

362.8651 

IV v,?t 

214.6-6? 
143. 0655 
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Figure A. 5. 29 
Model #2 
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Figure A. 5 . 31 
Model #? 
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Figure A. 5. 32 
Model ft 2 
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Figure A. 5. 33 
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Figure A. 5. 37 
Model #2 
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Figure A. 5 . 38 
Model #2 
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Figure A. 5.44 
Model #2 
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The size of the propellers were calculated using 
the following equation from Reference 24: 

1 _ n = + J{\ + ( T/A) - ( 2pV 2 ) } (A. 6.1) 


where: 


n 

T 

A 

P 

V 


propeller efficiency (85%) 
required Thrust 
actuator disc area 

density of air at cruising altitude 
cruising velocity 


With the values entered for thrust, density and 


velocity, a circular cross-sectional area A of the 
propeller disc as described in Reference 24. The 
effective area for each of the propellers is half of 
this and the resulting propeller diameter is 10* feet 
The rotational speed of the propeller was determined 


using the following equation: 


M* = J{ ( V/a) 2 + (nR/a)M (A. 6. 2) 

where : 

- maximum Mach number permitted by 
propeller tip = 0.8 

y _ forward velocity of aircraft 

n - rotational speed of propeller 

r - radius of propeller 

a - speed of sound 

Reference 11 described a procedure that produces an 
estimate of the size of a gearbox. Gearbox weight versus 
a factor identified as Q is plotted in Figure A. 6.1. Each 
plot is represented by another factor called K which is 
selected from the table in Figure A. 6. 2. 


A. 6. 2 



3 

Horsepower * (m + *) (A. 6. 3) 

Q ~ Pinion RPM m 


where : 

Pinion RPM - equivalent to n 

Horsepower - power produced by motor 

m - gear reduction ratio 

Once Q was calculated, K was selected as 600 

pounds per square inch (for epicyclic spur gears in 

aerospace aplications ) . These two values were cross- 

referenced to produce a weight of 80 pounds which was 

rounded up to 100 pounds since this method of sizing 

gearboxes is not very accurate. 


A. 6. 3 


Figure A. 6.1 
K-Factor Table 
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Figure A. 6. 2 

0-Factor vs. Gearbox Weight 



Gear weights plotted against Q factor for dtfferen* 
tensities of tooth loading (English uniu.i 
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Appendix A. 7 


Static Stability D erivation and Analysis Cod e 



Contribution of Aircraft Components 


It is of interest to know the contribution of the 
wing, fuselage, horizontal tail, and canard to the 
pitching moment and static stability characteristics of 
the airplane. Our emphasis will rely on methods that can 
be derived from simple theoretical considerations. These 
methods are generally accurate for the purposes of a 
preliminary design such as this. They show the 
relationship between the stability coefficients and the 
geometric and aerodynamic characteristics of the airplane. 
(Reference 23, p.44) 

Wing Contribution 


The contribution of the wing to the vehicle's static 

stability can be examined with the aid of Fig. A. 7.1. In 

this sketch the wing has been replaced by its mean chord. 

The distances from the wing leading edge to the 

aerodynamic center and the CG are denoted by X ac and X CG , 

respectively. The vertical displacement of the CG is 

denoted by The angle the wing chord line makes with 

J CG 

the fuselage reference line is denoted as i w * 


Summing moments about the CG, we obtain the following 
equation : 


A. 7. 2 


(A. 7.1) 


2 Moments = M CGw 

"cgw ’ Incest o ( ,-i w )tX CG -X ao ] + D u si n (a w -i tt ) 
t x CG‘ X ac ] + L w sin(0 M _i w UZ CG ] 

- D w cos(a w -i w )[Z CG ] + M ac ^ 

Dividing by ^pV 2 Sc yields: 

^CGW = C Lw (X CG /C - X ac/ C > 

(X CG /c " X ac /c) + C Lw Z CG /C 

sin(a w - i w ) - c jy v z cG^ C cos *“w - 

+ C Mac w 

Equation A. 7. 2 can be simplified by assuming that the angle 
of attack is small, with this assumption the following 
approximations can be made. 

cos(a w -i w )=l, sin(a w -i w )=a w -i w , C L »C D 

If we further assume that the vertical contribution is 
negligible, then Eq. A. 7. 2 reduces to: 

C M CG "’ =C Mac W+C Lw (X CG /C " X ac /C) (A. 7. 3) 


or 

c M CGw= c Mac « + ( c L Ow +dC L /dQ! w* a w ) ( X CG /C_X ac /C ) 


A. 7. 3 



where C LW -C L0w +dC L /da w a w . Applying the condition for 
static stability yields: 

C MO W=C Mac w+C LOw ( X CG /C-X ac /C * 

dC M /da w =dC L /da w ^ x cG^ c-X ac^ c 

Canard and Tail Contribution 

The horizontal tail and canard are analyzed in much the 
same way as the wing. The contributions of these components may 
be developed with the aid of Fig. A. 7. 2. Both the tail and 
canard were designed as all moving control surfaces to increase 
their effectiveness and decrease overall drag, (see Figure 
A. 7. 3) The angle of attack at the tail can be expressed as 
a _i w — £ +i t where e and i t are the downwash and tail 

incidences, respectively. The angle of attack at the canard may 
be written e c -o w -i w +i c where i c in the canard incidence angle. 
Note that the upwash at the canard has been assumed small and 

neglected . 

If we assume small angles and neglect the drag 
contribution of the tail, the total lift of the wing and 


(A. 7. 5) 
(A. 7. 6) 


tail can be expressed as: 


or 


C L C Lw 


+ ri s t c Lt /s 


(A. 7. 8) 


where 


^ = (Jspv t l )/(^pv w l ) ' * 

and ^ is the tail efficiency. (For this analysis, is 
assumed to equal unity). 


The pitching moment due to the tail can be 
obtained by summing the moments about the CG: 

M = -1 [L < .cos(arRL-6)+D t sin(aFRL-e) ] (A. 7. 

t t "C 

-Z„„t [ D+.COS ( aFRL- 6 ) -L t s in ( aFRL- e ) ] 

CG t. u 

Assuming again that C L »»C D and neglecting all but the 
first term, Eq. A. 7. 10 reduces to 


M t -_ 1 t L t _ “ 1 t C Lt^ pV t * S t 

C Mt =M t /!spVJSC = - 1 t S tf[ C Lt /SC 
C Mt = " V H^ C Lt 


A. 7. 5 



where v H -l t S t /Se, the horizontal tail volume ratio. The 
coefficient C Lt can be written as 

C Lt =(d V d 0 t>“t =(<, V d 0 t) <“w' i w‘ e+i t ) ( A. 7 . 14 ) 

where dC L /da t is the slope of the tail lift curve. The 
downwash angle e can be expressed as 

„ (A. 7. 15) 

e=e Q +( de/da) a w 

where £ Q is the downwash at zero angle of attack. 

The downwash behind an elliptically loaded wing can 
be derived from finite wing theory. 

(A. 7. 16) 

c=2C Lw /,rAR w 

where the downwash angle is in radians. The rate of 
change of downwash angle with angle of attack is 
determined by the derivative of Eg. A. 7. 16: 


( de /da ) =2 ( dC L /da w ) AAI^ 

These expressions do not take into account the relative 
positions of the tail plane relative to the wing. A more 
accurate method would require the fabrication of a model 
for testing in a wind tunnel. Rewriting the tail 


A. 7. 6 


contribution to the pitching moment yields 


(A. 7. 18) 

C^CGt=--V H f|C Lt 

w J _. • x (A. 7. 19) 

C M CGt=I l V H (dC L /da t ) ( e 0 +1 w" 1 t ) 

-^y R (dC L /da t )a (1-de/da) 

Finally, the expressions for the intercept and slope 
may be written: 


C MO t=f l V H (dC L /da t ) ( e O +i w _i t ) 

d c M /do t = -r^v H (dC L /de t )(l-de/de) ( A - 

By similar analysis, the intercept and slope of the 
canard are found to be: 

C M0 C= 1 V H (dC L /d< ‘c )(i w‘ i c> (A ‘ 

dC H /da c“1 V H (dC l. /d “c ) <A ' 

H.1.3: Fuselage contribution 


The fuselage contribution to the pitching moment 
curves is outlined in Ref. 23, pp. 49-51. The fuselage 
is divided into segments and the local induced angle due 
to the wing upwash or downwash for each segment can be 


A. 7. 7 



estimated. A computer program was written to estimate 
the intercept and slope of the fuselage contribution to 
the pitching moment curves; this code is included in 
this appendix. 

Stick Fixed Neutral Point 


The total pitching moment for the airplane can now 
be obtained by summing the wing, fuselage, and tail 
contributions . 

C H CG=C M0 +(dC M /da)£ ' (A. 7. 24) 

where 

C MO =C M° S+C M M+ (V V Ht (dC L/ d<I t> <£ 0 +1 w- i t ) (A. 7. 25) 

+ lV d V d<, c Hi w- i c> 

dC H /da=(dC L /da u ) ( X CG /c-X ac /c)+(dC H /do f ) ( A. 7 . 26 

-a V Ht (dC L /da t> ( 1 _de / da)+l \ V Hc (dC L /da c ) 

Setting dC M /da equal to zero and solving for the center 

of gravity position yields 

<x Hp /c)=( x ac /c ) - ( ( dC M /do £ ) / ( dC L /da„) )+([V H ( A . 7 . 27 

( (dC L /da t )/(dC L /da w ) ) (1-de/da)- 
f|V H ( ( dC L /da c ) / ( dC L /da w ) ) 


A. 7. 8 



This location, the stick fixed neutral point, is the 
point of neutral stability; movement of the CG behind 
this point will cause the airplane to become unstable. 

Static Stability Analysis Code 

The stability analysis of the aircraft required the 
creation of a computer program to quickly solve the 
longitudinal stability equations. This would allow us to 
see the effects of changes of variable values on the 
static stability of the aircraft, and thus to choose the 
best control surface configuration for the given flight 
conditions . 

The static stability analysis program performs four 
interrelated functions. The first is the calculation of 
the moment coefficients due to the wing, tail, and 
canard, as well as calculating the slope and intercept of 
Cm vs. a for each component. The second function is the 
optional calculation of the Cm, and slope and intercept 
of Cm vs. a for the fuselage of the aircraft. The third 
function is the tabulation (formatted for Lotus) and 
output of the results, for the components and the 
composite aircraft, vs. a. Finally, the program 
calculates the neutral point of the aircraft for 
subsequent iterations. 

All necessary geometric and performance constants are 
valued within the program. The distances to the center 


A. 7. 9 


of gravity and aerodynamic center (measured from the 
quarter chord of the wing), and tail and canard quarter 
chords (measured from the leading edge of the wing) , are 
left as variable inputs. Also left as variable inputs 
are the angles of incidence of the wing, tail, and 
canard. 


A. 7. 10 


3RCT** static stability Analysis ^ 

4 REM ** Program 

•si: .s!ss.s.s:.:-n25.s5.KK.s: :: 

8 REM 

20 R PRINT "This version, STAT7.BAS, accepts degrees as input." 

30 PRINT 
100 REM 
110 REM 
120 REM 
130 REM 
140 REM 

150 LET C =15.36 


Static Stability analysis program 


Geometric constants 


REM wing mean aerodynamic chord : 

Ho LET sw = 3594.25 : REM wing planform area ■ ^Sw 
170 LET ARW = 15.23 : REM wing aspect ratio . - ,ar 

REM tail planform area . • St 

REM tail aspect ratio : 

REM canard planform area 
REM canard aspect ratio . - ---- - = „. ifus 

• REM fuselage incidence angle . • i ' 

LET LC = 23.66: LET XAC=3.84: LET XCG-3.049 


?”ARt =”;art 
: ?»sc =” ;sc 
?”ARc =” ;arc 


? " ETA =" ;eta 


180 LET ST = 520 
190 LET ART =5.2 
200 LET SC = 470 : 

210 LET ARC =4.7 
220 LET IFUS = 0 
230 LET LT = 36.5 
240 REM 
250 REM 

260 REM performance constants 

7-n LET ETA = 1 ' REM wing/tail f if cf w^ 

III S lift-curve slope for infnite wing = T-Claa - 

320 CLAW = CLAA/ ( 1+ ( CLAA/ ( 3 . 14 159* ARW) ) ) 

330 CLAT - CLAA/ ( 1+ (CLAA/ ( 3 . 14159*ART) ) ) 

340 CLAC = CLAA/( l+(CLAA/( 3 . 14159*AR )^)) 

360 LET A0W = =*-1.5 fL^ot the wing at 0 lift : ? "A0w =";a0w 

370 REM 
380 REM 
390 REM 

400 REM Semivariable inputs 

420 INPUT "It ; " ; IT: IT=IT*1 . 745329E-02 
430 INPUT "Ic ;IC:IC=IC*1.745329E-02 
440 INPUT "Iw : ” ; IW: IW=IW*1 . 745329E-02 
450 REM 
460 REM 
470 REM 

480 REM , ^ . 

490 REM Preliminary calculations 

500 REM 

510 VHT = ( LT*ST ) /( SW*C) 

520 VHC = -(LC*SC)/(SW*C) 

530 DEDA = ( 2*CLAW) / ( 3 . 14 159* ARW) 

1 EO = (2*CL0W)/(3.14159*ARW) 

M REM 
560 REM 
570 GOTO 760 

580 REM Find CMaf and CMOf 


;claa 


CMOc , CMac 


600 INPUT "Do you need to calculate CMf »Q$ 

610 IF Q$="Y" THEN GOTO 650 
620 INPUT "CMOf :";CM0F 
63'' INPUT "CMaf : " ;CMAF 

tio INPUT 7 "How many sections?” .-SECTIONS 

660 CMAF=0 : CM0F=0 
/■on pot? t = 1 TO SECTIONS 

680 PRINT "Length of section" ;IsINPUT DXF 
(■on TNPUT "Wf for this section? ,wf 
700 INPUT "dEu/dA for this section?" ;DEUDA 
710 LET CMAF=CMAF+( (WF A 2 ) *DEUDA*DXF) 

720 LET CM0F=CM0F+(WF-2)*(A0W +I FUS)*DXF 

730 NEXT I 

740 CMAF = CMAF/ ( 36 • 5*SW*C ) 

750 CMOF — CM0F*DK/ ( 36 . 5*SW*C ) 

760 LET CM0F=-. 00158 
770 LET CMAF=. 04091 
780 REM 

800 REM Find CMOw , CMaw , CMOt , CMat 
Q 1 A REM 

820 CMOW = CMACW+CLOW* ( XCG/C-XAC/C ) 

830 CMAW = CLAW* (XCG/C-XAC/C) 

840 CMOT = ETA*VHT* ( CLAT* ( E0+IW-IT) ) 

850 CMOC = VHC* ( CLAC* ( IW-IC ) ) 

860 CMAT =-ETA*VHT* ( CLAT* ( 1-DEDA) ) 

870 CMAC =-VHC* ( CLAC ) 

880 REM 

9U REM Find CMw , CMf , ' T \ F n+TW-ir?+( ETA*VHC*CLAC) *( IW-IC ) 

HI I S:^;^^:^^(^T.ciAT,.,I-OED A )-(ETA.VHC.CLAC, 

940 REM 

960 A= .1745 

980 CMW = CM0W+CMAW*A 

990 CMF = CM0F+CMAF*A 

1000 CMT = CM0T+CMAT*A 

1010 CMC = CM0C+CMAC*A 

1020 CMCG = CM0+(CMA*A) 

1040 PRINT 

1070 PRINT" CM0=";CM0 

1080 PRINT" CMa=";CMA . 9 ,, reaion : " 7 ( CMA - CMAC) 

1085 PRINT "CMa in \tl™) region: "; (CMA - CMAT) 

1086 PRINT "CMa m non linear (tail) regxun 

1090 REM 

1100 REM Find Xnp 

XNP=C* ( XAC/C-CHAF/CLAW+ETA.VHT. ( CLAT/CLAW ) • ( 1-DEDA ) +VHC* CLAC /CLAW ) 

1133 PRINT"Stability Margin = " : ( (CMCG-CM0)/1 .184) 

1140 REM 

1150 REM , . curves at 0 AOA and 10 AO A" 

1160 PRINT "Sample of points on curve 10 DEG REES" 

11-70 PRINT "COMPONENT: 0 DEGREES 


1 3 PRINT 

li90 PRINT "Wing 
1200 PRINT "Canard 
1210 PRINT "Tail 


" ; CMOW ; TAB (40) CMW 
" ; CMOC ; TAB ( 40 ) CMC 
" ; CMOT ; TAB (40) CMT 



1220 

1230 

1240 

1250 


PRINT "Fuselage " ;CMOF;TAB( 40)CMF 

PRINT "Airplane " ;CM0 ;TAB( 40 )CMCG 

PRINT 
RUN 
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Figure A. 7. 2 

Horizontal Tail Contribution to the 
Pitching Moment 
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Appendix A. 8 
Dynamic Stability, 



The example on p. 131 of Reference 23 was used as a 
test case for the block diagram input into ASDEQ to 
solve equation 7.2.1. The theoretical approximations of 
the problem were compared to the results of the program 
(modified to show the long and short periods separately 
neglecting interactions). This comparison is summarized 
in the following table. 

* 

Theoretical approx imations Results of proqra jn 

Psp=2 . 26 sec Psp=2 . 3 sec 

Pp=25 . 54 sec Pp=24.5 sec 

* Extrapolated from Figure A. 8.1 & A. 8. 2 

The unmodified program, which takes the 
interconnectivity of the modes into account, was tested 
by solving the following equations simultaneously using 
data from Figure A. 8. 3. The graphical results are 
presented in Figures A. 8. 4 & A. 8. 5. 

p p=-(l/T d )ln(u 2 /u 1 ) 
p(l- p 2 ) 1/2 =2 /T d 

The results compare to the example's exact results 
as follows. 

Rpferenr^ results pypqram result s 


A. 8. 2 



p= . 2 1 4 4 rad/sec 


p= . 21 rad/sec 


The above comparisons validate the accuracy of both 
the main block diagram and the long and short period 
modifications . 


A. 8. 3 



Figure A. 8.1 
Short Period Mode 




Figure A. 8. 2 
Long Period Mode 
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Figure A. 8. 3 
Combined Mode Data 



Data for this run was read from ft : DVNBOOK. PSD 
under automatic step control with local error — . 
or. 04- 17- 1989 at . £0:50:45 

IND. VfiR BLOCK 1 BLOCK 8 BLOCK 18 

du/dt dw/dt dq/dt 


0. 000 

0 

0. 500 

3. 888932E— 04 

1 . 000 

1. 1 31991 E— 03 

1 . 500 

1. 958666E-03 

2. 000 

2. 73988E— 03 

2. 500 

3. 467492E— 03 

3. 000 

4. 143304E— 03 

3. 500 

4. 762 1 58E-03 

4. 000 

5. 31 7657E— 03 

4. 500 

5. 804466E— 03 

5. 000 

6. 21 8296E-03 

5. 500 

6. 555759E— 03 

6. 000 

6. 814365E— 03 

6. 500 

6. 992543E— 03 

7. 000 

7. 089645E— 03 

7. 500 

7. 10594 IE— 03 

8. 000 

. 0070426 

8. 500 

6. 90 166E-03 

9. 000 

6. 68599 IE-03 

9. 500 

6. 39925E-03 

10. 000 

6. 045828E— 03 

10. 500 

5. 630785E— 03 

1 1 . OOO 

5. 159786E— 03 

1 1 . 500 

4. 639028E— 03 

12. 000 

4. 075157E— 03 

12. 500 

3. 475193E— 03 

13. 000 

2. 846438E— 03 

13. 500 

2. 196396E-03 

14. 000 

1 . 532683E— 03 

14. 500 

8. 629402E— 04 

15. OOO 

1. 947489E— 04 

15. 500 

-4. 64454 1 E— 04 

16. OOO 

-1. 107457E— 03 

16. 500 

-1. 727347E-03 

17. OOO 

-2. 31 7589E-03 

17. 500 

-2. 872085E-03 

18. OOO 

—3. 385239E— 03 

18. 500 

-3. 85201 IE— 03 

19. OOO 

-4. 267962E— 03 

19. 500 

-4. 629295E— 03 

20. OOO 

-4. 932888E— 03 

20. 500 

-5. 176314E— 03 

21. OOO 

-5. 357859E— 03 

21. 500 

-5. 476531E— 03 

22. OOO 

—5. 532053E— 03 

£2. 500 

-5. 52486 IE— 03 

23. OOO 

-5. 456083E— 03 

23. 500 

-5. 32751 5E-03 

£4. OOO 

-5. 141594E-03 

£4. 500 

-4. 901353E— 03 

25. OOO 

-4. 6 1 0383E— 03 

25. 500 

-4. 272782E— 03 

£6. OOO 

-3. 693099E— 03 

26. 500 

—3. 476278E— 03 


.017 

. 00£1783 
-1. 008656E-03 
-4. 4926£3E— 04 
-1. 519205E-04 
-1. 82749E-04 
-2. 458263E— 04 
-2. 856293E-04 
-3. 159344E— 04 
-3. 426886E— 04 
-3. 656982E— 04 
-3. 842862E— 04 
-3. 9824 1 8E— 04 
-4. 075071E— 04 
-4. 12067 IE— 04 
-4. 1 19495E— 04 
-4. 072329E-04 
-3. 980467E— 04 
-3. 845682E— 04 
-3. 670198E— 04 
-3. 456654E— 04 
-3. £0807 IE— 04 
-2. 9278 IE— 04 
-2. 61953E-04 
-2. 287144E-04 
-1. 934767E-04 
-1. 566668E-04 
-1. 18722E— 04 
-8. 008524E— 05 
-4. 1 19949E-05 
—2. 503266E— 06 
3. 557443E— 05 
7. 261863E-05 
1. 082329E-04 
1. 420433E— 04 
1. 737024E-04 
2. 028925E-04 
2. 293288E— 04 
2. 527622E-04 
2. 729809E— 04 
2. 8981 31 E— 04 
3. 03 127 IE— 04 
3. 128331E-04 
3. 1 88829E— 04 
3. 2 1270 IE— 04 
3. 200293E— 04 
3. 15234BE-04 
3. 069996E-04 
2. 95473 IE-04 
2. 80839 IE-04 
2. 63313E-04 
2. 43 1 387E-04 
2. 2O506E— 04 
1 . 959464E— 04 


0 

6. 7835E-05 
1. 044966E-05 
5. 3870B6E— 06 
4. 652566E— 06 
4. 277175E— 06 
5. 170457E— 06 
6. 179634E— 06 
7. 034185E— 06 
7. 769967E— 06 
8. 4076E-06 
8. 942817E— 06 
9. 369498E— 06 
9. 684288E— 06 
9. 885614E— 06 
9. 97317E— 06 
9. 947877E— 06 
9. 81 1899E-06 
9. 568606E-06 
9. 222508E— 06 
8. 77918BE— 06 
8. 24522E— 06 
7. 628083E— 06 
6. 936061E— 06 
6. 178132E— 06 
5. 363862E— 06 
4. 503289E— 06 
3. 606797E— 06 
2. 68499BE-06 
1. 74B612E-06 
8. 083396E— 07 
-1. 252517E— 07 
-1. 041847E-06 
-1. 931492E-06 
-2. 784702E— 06 
-3. 592555E-06 
-4. 346787E-06 
-5. 039868E-06 
-5. 665082E— 06 
-6. 216578E— 06 
-6. 689427E— 06 
-7. 079664E— 06 
-7. 3843 IE— 06 
-7. 60 139 IE— 06 
-7. 729946E-06 
-7. 770021E— 06 
-7. 722642E— 06 
-7. 589803E— 06 
-7. 374414E— 06 
-7. 0B0258E— 06 
-6. 71 194 IE— 06 
-6. £748 IE— 06 
-5. 77490 IE-06 
-5. 218845E— 06 



£8. OOO - 
£8. 500 - 

£9. 000 
£9. 500 ■ 

30. 000 

30. 500 

3 1 . 000 

31. 500 
3£. OOO 
3£. 500 

33. 000 
33. 500 

34. 000 

34. 500 

35. 000 

35. 500 

36. 000 

36. 500 

37. 000 

37. 500 

38. 000 

38. 500 

39. 000 

39. 500 

40. 000 

40. 500 

4 1 . 000 

41. 500 

42. 000 

42. 500 

43. 000 

43. 500 

44. 000 

44. 500 

45. 000 

45. 500 

46. 000 

46. 500 

47. 000 

47. 500 

48. 000 

48. 500 

49. 000 

49. 500 

50. 000 

50. 500 

51 . OOO 

5 1 . 500 

52. OOO 

52. 500 

53. OOO 

53. 500 

54. 000 

54. 500 

55. 000 

55. 500 

56. 000 

56. 500 

57. 000 

57. 500 

58. 000 

58. 500 

59. 000 


,'=i. sue > - 


£. 05701 86— 03 
1 . 546747E— 03 
1 . 027726E-03 
5. 05896 IE -04 
1. 28621 4E-05 
5. 2281 19E-04 
1 . 018403E— 03 
1. 494339E— 03 
1. 945634E-03 
2. 367659E-03 
2. 75619E-03 
3. 1 07445E— 03 
3. 418126E— 03 
3. 685439E— 03 
3. 907 1 25E-03 
4. 081468E-03 
4. 207314E— 03 
4. 284067E-03 
4. 31 1689E-03 
4. 290693E-03 
4. 222125E— 03 
4. 1 07547E-03 
3. 949007E— 03 
3. 749009E— 03 
3. 5 1048 IE— 03 
3. 236729E-03 
2. 93 140 IE— 03 
2. 598434E— 03 
2. 242008E— 03 
1. 866495E— 03 
1. 47641 IE-03 
1. 076353E— 03 
6. 709574E— 04 
£. 648409E— 04 
-1 . 374485E— 04 
-5. 314818E— 04 
-9. 12996E-04 
-1. 277939E-03 
-1. 62251 2E-03 
-1. 943207E-03 
-£. 23684E— 03 
— 1~* _ 500585E— 03 
-2. 731995E— 03 
-2. 929025E— 03 
-3. 09005E— 03 
-3. 213872E— 03 
-3. 29973 IE— 03 
0033473 
-3. 35669E-03 
-3. 328432E— 03 
—3. 263472E— 03 
-3. 163148E-03 
-3. 0291 7 IE— 03 
0028636 
-2. 6688 1 2E-03 
-£. 44747 1 E— 03 
-2. 202493E— 03 
-1. 937008E— 03 
— 1 . 654324E— 03 
-1. 357886E-03 
- 1 . 05 1 235E-03 
-7. 379656E— 04 
-4. 216879E-04 

1 |-)S‘=»A4AP-04 


1. i£t.76t.-04 
6. £91 658E-05 
5. £7£S£9E-05 
£. £455422—05 
-7. 5624Q6E— ‘>6 
—3. 699442E-05 
-6. 552146E-05 
-9. 284009E— 05 
-1. 1 86654E-04 
-1. 427339E-04 
-1 . 64806 IE-04 
-1. 846688E-04 
-£. 021372E-04 
-£. 17056 IE- 04 
-£. 2930 1 5E— 04 
-2. 387812E-04 
-£. 454357E— 04 
-£. 49237 BE— 04 
-2. 501926E-04 
-2. 483372E-04 
-£. 43739 IE-04 
-£. 364955E— 04 
— £. 267316E— 04 
-2. 145984E— 04 
-£. 00271 IE-04 
-1 . 839465E-04 
-1. 6584E-04 
-1. 461839E— 04 
-1. £522 34 E -04 
-1 . 032145E— 04 
-8. 042053E— 05 
-5. 7 1 09 1 8E-05 
-3. 354944E— 05 
- 1 . 000859E— 05 
1 . 325074E— 05 
3. 597348E— 05 
5. 791481E-05 
7. B84274E— 05 
9. 854043E— 05 
1. 168085E— 04 
1 . 334667E— 04 
1. 483559E— 04 
1. 6 1 3394E— 04 
1 . 72304E— 04 
1. 81 161 E— 04 
1 . 878467E— 04 
1. 92 3227 E— 04 
1. 945758E-04 
1. 946161E-04 
1. 92486 IE-04 
1. 882396E— 04 
1 . 81962E-04 
1. 737568E-04 
1 . 637481E— 04 
1 . 520777E— 04 
1. 389037E— 04 
1 . 243984E— 04 
1. 08746 IE-04 
9. 214067E— 05 
7. 478359E— 05 
5. 688 1 £7 E— 05 
3. 864272E— 05 
£ . 02773E-05 
1 _ -Of-. 


-3. 2873288 -08 
-£. 581986E- 08 
-1. 85959 IE- 06 
-1. 128495E- 06 
-3. 970101E— 07 
3. 2668B6E-07 
1 . 0346462-06 
1. 719223E— 06 
£. 37316E-06 
£. 969666E-06 
3. 56248E-06 
4. 085936E— 06 
4. 55501 IE -06 
4. 965376E- 06 
5. 31 3427E— 06 
5. 596318E 06 
5. 81 1979E-06 
5. 959124E— 06 
6. 037254E— 06 
6. 046649E— 06 
5. 988353E— 06 
5. 864149E- 06 
5. 676523E-06 
5. 428632E-06 
5. 12425E-06 
4. 767723E— 06 
4. 36390 IE -06 
3. 918067E —06 
3. 43596E-06 
£. 923512E— 06 
£. 386974E— 06 
1. B32739E-06 
1. £67£95E-06 
6. 971433E-07 
1. 287307E-07 
-4. 316231E— 07 
-9. 777946E— 07 
-1. 50392E-06 
-£. 004459E-06 
— £. 47424BE- 06 
-£. 908553E— 06 
—3. 3031 1 BE-OE. 
-3. 654198E-06 
-3. 9586E-06 
-4. £1 3703E-06 
-4. 417482E— 06 
-4. 568519E— 06 
-4. 666007 E~ 06 
-4. 709752E— 06 
-4. 700164E- 06 
-4. 63B241E-06 
-4. 525548E— 06 
-4. 364193E— 06 
-4. 156791E— 06 
-3. 906428E— 06 
—3. 61661 7E— 06 
3. 291254E-06 
£. 934567E-06 
£. 551064E— 06 
£. 145476E— 06 
1. 72270 IE- 06 
•1. 287749E- 06 
■8. 4568 IE— 07 
■ h . i'll SS34P- i'. 7 





tL l . OoO 

£1 . 500 
62. 000 
62. 500 
£3. 000 
£3. 500 
£4. 000 
£4. 500 
£5,. 000 
£5. 500 
££. 000 
66. 500 

67. 000 

67. 500 
68. 000 

68. 500 

69. 000 

69. 500 
70. 000 

70. 500 

7 1 . 000 

7 1 . 500 

72. 000 

72. 500 

73. 000 
73. 500 

74. 000 

74. 500 
75. 000 

75. 500 

76. 000 

76. 500 

77. 000 

77. 500 

78. 000 

78. 500 

79. 000 

79. 500 

80. 000 

80. 500 
8 1 . 000 

8 1 . 500 
82. 000 

82. 500 

83. 000 

83. 500 

84. 000 

84. 500 

85. 000 

85. 500 

86. 000 

86. 500 

87. 000 

87. 500 

88. 000 

88. 500 

89. 000 

89. 500 

90. 000 

90. 500 
91. OOO 

9 1 . 500 

92. 000 

92. 500 


ti . -V-J i H l- " 1 H- 

1 . 0S25++E-03 
1 . 345246E-03 
1 . 588538E-03 
1 . 81 0024E— 03 
2. 007593E— 03 
2. 179439E— 03 
2. 324072E— 03 
2. 440332E— 03 
2. 527393E— 03 
2. 584772E-03 
2. 612324E— 03 
2. 610238E-03 
2. 579034E— 03 
2. 519547E— 03 
2. 432913E— 03 
2. 320553E— 03 
2. 184153E— 03 
2. 025637E 
1. 847145E-03 
1. 651005E— 03 
1. 439702E— 03 
1. 215849E— 03 
9. 82157E— 04 
7. 41 3993E— 04 
4. 963834E— 04 
2. 49917E-04 
4. 776952E— 06 
■2. 363209E-04 
■4. 70751 3E-04 
-6. 960074E— 04 
■9. 097265E— 04 
-1. 109715E— 03 
■1. 293967E— 03 
■ 1 . 460689E— 03 
■1 . 608312E— 03 
-1. 735505E-03 
-1. 841 189E-03 
-1. 924544E-03 
■1. 985012E-03 
-2. 022299E-03 
-2. 036376E— 03 
-2. 027474E— 03 
-1. 996076E-03 
-1. 942909E-03 
-1. 86893E-03 
■1. 77531 4E-03 
-1 . 663433E-03 
-1. 534843E-03 
-1. 391258E— 03 
-1. 234534E-03 
-1. 06664E-03 
-8. 896376E— 04 
-7. 056562E— 04 
-5. 16867 IE— 04 
-3. 254594E— 04 
-1. 336145E— 04 
5. 651762E-05 
2. 42842E-04 
4. 23341 IE— 04 
5. 960965E— 04 
7. 593079E— 04 
9. 1 13121E-04 
1 . 050598E— 03 


-£. £ 39754 E- 05 
-8. 13991 3E— 05 
-9. 524095E-05 
-1 . 077878E— 04 
-1. 18921 IE-04 
-1 . 285403E-04 
— 1 . 36563E— 04 
-1. 429265E— 04 
-1. 475874E-04 
— 1. 50522E— 04 
-1. 517264E-04 
-1. 512161E-04 
-1. 490254E— 04 
-1. 452069E— 04 
-1 . 398303E— 04 
-1. 32981 9E-04 
-1. 247626E-04 
— 1. 152872E— 04 
-1. 046824E— 04 
-9. 308567E— 05 
-8. 0643E-05 
-6. 75075 IE-05 
-5. 383751E— 05 
-3. 979458E— 05 
-2. 554181E— 05 
-1. 124188E— 05 
2. 94469E— 06 
1 . 686 1 36E— 05 
3. 035721E-05 
4. 328852E— 05 
5. 55203 IE— 05 
6. 692766E— 05 
7. 73969E— 05 
8. 66ii'66E 05 
9. 512943E— 05 
1 . 022308E— 04 
1 . 0807 1 8E-04 
1. 12608 IE— 04 
1. 1581 IE-04 
1. 176669E— 04 
1. 181775E-04 
1. 173598E— 04 
1 . 152449E— 04 
1. 1 18782E-04 
1. 073182E-O4 
1. 016357E-04 
9. 491306E— 05 
8. 724264E— 05 
7. 872596E— 05 
6. 947233E— 05 
5. 959733E— 05 
4. 922169E-05 
3. 846965E— 05 
2. 746759E— 05 
1. 634254E— 05 
5. 22084 IE-06 
-5. 773366E— 06 
-1.651941 E— 05 
-2. 690142E— 05 
—3. 68095 IE— 05 
—4. 614096E— 05 
-5. 4801 16E-05 
-6. 270462E-05 


1 . 29/71 4t- c 
1 . 6603 1 4E— Ofc 
2. 037754E- 06 
2. 366464t-<.»£ 

2. £6325 IE— 06 
2. 925334E— 06 
3. 150366E— o£ 

3. 336453E— 06 
3. 482166E— o£ 
3. 586552E— 06 
3. 649138E-06 
3. 66992 IE-06 
3. 649368E— 06 
3. 5884E-06 
3. 48837 IE— 06 
3. 351053E— 06 
3. 1 78601 E— 06 
2. 97 35 29 E— 06 
2. 736672E— 06 
2. 477152E— 06 
2. 192334E-06 
1. 887788E-06 
1. 567246E— 06 
1 . 234554E— 06 
8. 93633E-07 
5. 484263E— 07 
2. 0286 1 BE— 07 

-1. 39 195 IE— 07 
-4. 739827E— 07 
-7. 978B33E— 07 
-1. 10746 IE— 06 
-1. 399497E— 06 
-1. 671024E-06 
-1. 9 1935 IE— 06 
-2. 142092E— 06 
—2. 337188E— 06 
-2. 50292E— 06 
-2. 637926E-06 
-2. 741218E— 06 
-2. 812167E— 06 
-2. 85052 IE-06 
-2. 856398E— 06 
-2. 830273E— 06 
-2. 772975E— 06 
-2. 685662E— 06 
-2. 5698 1 3E— 06 
-2. 427197E— 06 
-2. 259S56E— 06 
-2. 070069E— 06 
-1. 860333E-06 
-1. 633324E-06 
—1. 391865E— 06 
-1. 138894E— 06 
-8. 774296E— 07 
-6. 105334E— 07 
—3. 41 2773E-07 
-7. 270802E— 06 
1. 921864E— 07 
4. 505087E— 07 
6. 994837E-07 
9. 364866E-07 
1. 159o7E— u6 
1. 36499E— 06 


omcmmi page is 
ot- poor quality 



93. 000 

1. 175824E— 03 

-6. 977568E— 05 

1 ■ 552224E— 06 

93. 500 

1 . 285824E-03 

-7. 594932E— 05 

1. 7 1 8994E— 06 

94. 000 

1. 379626E— 03 

-8. 1 171 SEE— 05 

1. 663777E— 06 

94. 500 

1. 456451 E -03 

-8. 539992E-05 

1. 985326E-06 

95. 000 

1. 515727E-03 

-8. 860382E— 05 

2. 08267 IE— 06 

95. 500 

1 . 557083E— 03 

-9. 07646E— 05 

2. 155128E— 06 

96. 000 

1. 580357E-03 

— 9. 1 87555E— 05 

2. 202304E— 06 

96. 500 

1. 58559 IE-03 

-9. 1 94 1 79E-05 

2. 224094E— 06 

97. 000 

1 . 573029E— 03 

-9. 098003E— 05 

2. 22068 IE— 06 

97. 500 

1. 543108E-03 

-8. 90 1 82 1 E— 05 

2. 1 92521 E— 06 

98. 000 

1. 496452E-03 

-8. 609494E— 05 

2. 140345E— 06 

98. 500 

1. 433862E-03 

8* 22589.1E — 05 

2. 065134E— 06 

99. 000 

1 . 356303E— 03 

-7. 756826E— 05 

1. 9681 16E-06 

99. 500 

1. 264892E-03 

-7. 20896E— 05 

1. 850736E-06 

100. 000 

1. 16088E-03 

- 6. 5897 SEE— 05 

1. 714647E-06 


Exact solution of example problem 


proof of program function 



dw/dt (V-axis) vs IHDEP.UAR (X-axis) 


Figure A. 8. 4 
Combined Mode 



Figure A. 8. 5 
Combined Mode 
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9.2 Pover Analysis 


9.2.1 Level Flight 



Power Analysis 


100 REM 
110 REM 

l£0 REM Level Flight 

1 REM 

1 REM ******* Variables and Constants ******* 


150 REM 
160 REM 
170 REM 
ISO REM 
190 REM 
£00 REM 
£ lO REM 
££0 REM 
£30 REM 
£40 REM 
£50 REM 
£60 REM 
£70 REM 
£80 REM 
£90 REM 
300 REM 
310 REM 
3i=! O REM 
330 REM 
340 REM 
350 REM 
360 REM 
370 REM 
380 REM 
3 REM 
4c«J REM 
410 REM 
420 REM 
430 REM 
440 REM 
450 REM 
460 REM 
470 
480 
490 
500 
510 

520 

530 

540 

550 

560 

570 

580 REM 
590 REM 
600 REM 
610 
620 
6’h 
6 

650 

660 

670 

680 


M 

S 

RO 

V 

AR 

FE 

NF 

L 

W 

CD 

CO 

K 

CL 

PI 

P 

D 

E 

PA 

RG 

EM 

EP 

PD 

FA 

PR 

T 

FR 

AI 

AD 


- Mach Number 

- Wing Area <ft' s 2) 

- Density (lb/ft^S) 

- Velocity (ft/s) 

- Asoect Ratio 

- Bank Angle (degrees) 

- Bank angle (radians) 

- Lift (lb) 

- Weight (lb) 

- Coefficient of Drag 

- Coefficient of Drag at zero lift 

- Constant 

- Coefficient of Lift 


- 3. 14 

- Power (hp) 

- Drag 

- Airolane efficiency 

- Power available (can be generated by the 

- Power generated /power received by motor 

- Motor efficiency 

- Prooeller efficiency 

- Power Density (hp/ft y '2) 

- Alt itude (ft ) 

- Power required (must be supplied to the 

- Thrust 

- Flight radius 

- Angle of incidence 

- Angle of incidence (degrees) 


motor 


mot or ) 


) 


******* Define Constants ******* 


E = 
PI 
I « 
ST 
K : 

G « 
RG 
EM 
AP 
CP 
PD 


. 85 

3. 1416 
. 01 
5 

.007 


32. 174 
.97 
. 98 
> 175 
550 


: REM Airplane efficiency 

• REM Value of one radian 
: REM Increase the Mach # 

• REM Increase the bank angle 

s rem Viscous Drag Coefficient estimated 
from a NACA 2214 airfoil 
: REM Gravitational constant (ft/s 2) 
s R£M Power recieved/ power generated 

• REM Motor efficiency 

: REM Prooeller area (ft' s 2) 

:REM Conversion from hp to ft*lbf/s 


<700* 1 . 3405* 1 0 ''-3* < . 3048 ) "'2 ) 


******* Data Entry ******* 

INPUT “Gross Weight" ;W 

INPUT "Wing Area"; S 

INPUT “Asoect Ratio" ;AR 

INPUT "Altitude of flight";FA 

INPUT "Soeed of sound at altitude ;SS 

INPUT "Density at altitude" ;R0 

INPUT "Coefficient of Drag © zero lift ;CO 

INPUT "Total propeller area";AP 


690 REM 

700 REM ******* Print Input Constants ******* 



’£0 

'30 

’40 

'50 

’60 

'70 

’80 

(00 REM 

(10 REM ******* 

(20 REM 

(30 

(40 REM 

(50 

(60 

(70 REM 
(80 

(90 REM 
(00 

(10 


(£0 

930 REM 
940 
950 
960 

970 REM 
960 

990 REM 

1 000 

L 0 REM 

1020 

1030 

1040 REM 

1050 

1060 

1070 REM 
1060 

1090 REM 
1 1 00 

1110 REM 
1120 


PRINT "Gross Weight „ 

PRINT "Wing Area .. ’qo 

PRINT "Aspect Ratio ? ° R 

PRINT "Altitude of flight 

PRINT "Speed of sound @ altitude " ;SS 

PRINT "Density at altitude 

PRINT "Coefficient of Drag 6 * zero lift 1 * 

PRINT "Total propeller area ”;AP 

Calculations ******* 

OPEN " A : 100CR2. out " FOR OUTPUT AS 1 LEN=£000 


L = W 

FOR M = .01 TO 1.1 STEP I 


V = M*SS 

CL = <L/(.5*R0*S*(V--'2> 

CD = (CO + K*(CL''2> + 

D = CD*. 5*R0* (V''2> *S 

P = D*V/CP 

PR = PD*S 

PA = RG*EM*EP*PR 

T = ( <PR/V)*CP) 


: REM VELOCITY 

) ) s REM COEFF OF 
LIFT 

(CL^'E) / <E*PI*AR) ) 

: REM COEFF OF 
DRAG 

sREM DRAG 

: REM POWER REQUIRED (ho) 
: REM POWER RECEIVED 
:REM POWER AVAILABLE 

: REM THRUST 


EP = (£/ ( 1+ ( 1+ (T/ (. 5*R0*(V) ''2*AP> ) 5) ) 

! REM PROP EFFICIENCY 


PA = RG*EM*EP*PR 
DH » < (PA - P) /W> *CP 


: REM POWER AVAILABLE 
iREM RATE OF CLIMB 


PRINT "MACH #";M;“ CLIMB RATE ";DH 
PRINT#1, M ; DH 


NEXT M 
CLOSE 1 
END 


A. 9. 2 
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9.2.2 Turning Flight 


Flight Preforrnance 
Turning Flight 


•a***#**********-* 


100 

110 

l £0 

130 


REM 

REM 

REM 

REM 

REM 


**************** 


Variables arid Constants 


160 

REM 

M 

- Mach Number 

170 

REM 

S 

- Wing Area (ft''£) 

ISO 

REM 

V 

- Velocity <ft/s> 

190 

REM 

AR 

- Aspect Ratio 

£00 

REM 

FE 

— Bank Angle (degrees) 

£10 

REM 

NF 

- Bank angle (radians) 

££0 

REM 

L 

- Lift (lb) 

£30 

REM 

W 

- Weight (lb) 

£40 

REM 

CD 

— Coefficient of Drag 

£50 

REM 

CO 

— Coefficient of Drag at zero lift 

£60 

REM 

K 

- Constant 

£70 

REM 

CL 

- Coefficient of Lift 

£60 

REM 

PI 

- 3. 14 

£90 

REM 

P 

— Power (hp) 

300 

REM 

D 

- Drag 

310 

REM 

E 

- Airplane efficiency 

3£0 

REM 

PA 

- Power available (can be generated by the motor) 

330 

REM 

RG 

- Power generated/Dower received by motor 

340 

REM 

EM 

- Motor efficiency 

350 

REM 

EP 

- Propeller efficiency 

360 

REM 

PD 

- Power Density (hp/ft"'£) 

370 

REM 

FA 

- Altitude (ft) 

380 

REM 

PR 

— Power reouired (must be supplied to the motor) 

3°0 

REM 

T 

- Thrust 

A 

REM 

FR 

- Flight radius 

410 

REM 

AI 

- Angle of incidence 

420 

REM 

AD 

- Angle of incidence (degrees) 

430 

440 

REM 

REM 

**************** Define Constant 5 ************************ 

450 

REM 


= .65 : REM Airolane efficiency 

460 


E 

470 


PI 

= 3.1416 : REM Value of one radian 

460 


I 

= .01 : REM Increase the Mach 4 

490 


ST 

= £ : REM Increase the bank angle 

500 


K 

= .007 : REM Viscous Drag Coefficient estimated 

510 


G 

from a NACA ££14 airfoil 

= 32. 174 :REM Gravitat ional constant (ft/s''£) 

5£0 


RG 

= .97 : REM Power recieved/ power generated 

530 


EM 

= .96 : REM Motor efficiency 

540 


AP 

s 175 : REM Prooeller area (ft A £) 

550 


CP 

= 550 : REM Conversion from hp to ft*lbf/s 

560 

570 

REM 

PD 

= (700*1. 3405*10'''— 3* (. 3048) ^2) 


560 

590 


REM 

REM 


**************** 


Data Entry 


******************************** 


600 REM 
610 
6£0 
630 
f 

6^rO 

660 

670 

680 

690 REM 


INPUT "Gross Weight" ;W 

INPUT "Wing Area" : S 

INPUT "Ascect Ratio" ;AR 

INPUT "Altitude of flight" ;FA 

INPUT "Soeed of sound at alt itude" ;SS 

INPUT "Density at altitude" $R0 

INPUT "Coefficient of Drag 0 zero lift";C0 

INPUT "Total propeller area";AP A. 9. 9 


700 REM **************** Print Incut Constants' ************************ 


7£0 

LPRINT 

"Gross Weight 

** ; W 

730 

lprint 

"Wing Prea 

? s 

740 

lprint 

"Pspect Ratio 

" ;PR 

750 

lprint 

"Pltitude of flight 

" ;FP 

760 

LPRINT 

"Soeed of sound ® altitude 

" ;SS 

770 

lprint 

"Density at altitude 

" ; RD 

7°0 

LPRINT 

"Coefficient of Drag & zero lift 

" ;C0 

I 

LPRINT 

"Total propeller area 

" ;PP 


800 

810 

820 

830 


REM 

REM 

REM 


****************** 


Ca 1 cu 1 at i ons ****************************** 


OPEN "c:50CR. out 11 FOR OUTPUT PS 1 LEN=2000 


840 REM 
850 

860 REM 
870 
880 

890 REM 

900 

910 

920 

930 

940 

950 REM 
960 

970 REM 
980 

990 REM 
1 000 
1010 
1020 

1030 REM 
VO 
IvaO 
1060 


FOR FE = 3 TO 13 STEP ST 

NF = FE*PI / 180 
L = ( W/COS ( NF > > 

PRINT "Bank Pngle " ;FE; "degrees" 

PR I NT# 1, "Bank Pngle";FE 
PRINT 

PRINT "Mach #", "Vel ", "CL", "Cl imb Rate" 

PRINT 

FOR M = .01 TO 1.1 STEP I 
V = M*SS 

CL = (L/ (. 5*R0*S* (V''2> ) > 

CD = (CO + K* (CL^'E) 4- (CL''2> / <E*PI*PR> ) 
D * CD*. 5*R0* <V~2> *S 

FR = < V'2/ (G*TPN (NF) ) ) 

fil = PTN < V' s £/ (FR*G> ) + PTN (FR/FP) 

PD ■ PI* ( 180/ PI ) 


1070 
1080 
1090 
1100 
1110 
1120 
1130 
1 140 
1150 
1160 
1170 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 


REM 


P = D*V/CP : REM 

Power 

Required 



PR = PD*S*COS (PI ) : REM 

Power 

Received 

REM 


T = ( (PR/V)*CP) 



REM 


EP = (£/<l+(l + (T/(.5*R0*(V>' x £*PP) 5) ) 

REM 


PA = RG*EM*EP*PR sREM 

Power 

Pvai ladle 

REM 


DH * ( (PP - P) /W)*CP 



REM 


PRINT M, V, CL, PP, PR, DH 
PRINT#1, M : V ; CL ; PP ; P ; DH 



REM 

NEXT M 




REM 

NEXT FE 




REM 

CLOSE 1 




REM 

END 





(hp) 
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n> rj* n> t# 

(V li> Cv rji U 


Optimal Climb Conditions 

Chosen from Climb Rate Analysis 
Climb Rate maximized 


Wi «e M , velocity Clift 8? 'T 

“Trtr*‘‘“T 


25 

50 

75 

100 


0.09 

0.14 

0.22 

0.38 


91.44 

135.52 

212.97 

381.216 


. 0.416338 11.88128 2970.404 
5 0.555369 10.6274 6524.540 
7 0.744382 7.575022 11481.22 
9 0.802255 0.108634 28518.36 
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Hunter of turns 

(change in altitude per turn) 


dv/dh(xlO A 3) 1.2807 
dr/dh(xl0 A 3) 106.0047 
((dh/dt)/dh)*10 A 3)0.013974 


iltitude 
10 A -3 ft 
********* 

5 

6.023867 

7.157269 

8.407399 

9.781111 

11.28478 

12.92421 

14.70449 

16.62995 

18.70410 

20.92958 

23.30820 


Velocity Radius Onega Tine of Clinb change in 

ft/s ft Period ft/s altitude 

****************************************************** 
65.826 850.3101 0.077414 88.25075 11.6018 1.023867 
67.13726 958.8449 0.070018 97.57155 11.61610 1.133401 
68.58881 1078.990 0.063567 107.4738 11.63194 1.250130 
70.18985 1211.510 0.057935 117.9210 11.64941 1.373711 
71.94916 1357.130 0.053015 128.8648 11.66861 1.503673 
73.87492 1516.526 0.048713 140.2463 11.68962 1.639427 
75.97453 1690.313 0.044946 151.9979 11.71253 1.780281 
78.25454 1879.032 0.041646 164.0450 11.73741 1.925464 
80 72048 2083.140 0.038749 176.3085 11.76431 2.074149 
83.37684 2303.010 0.036203 188.7073 11.79330 2.225482 
86.22702 2538.921 0.033962 201.1613 11.82440 2.378611 
89 27331 2791.065 0.031985 213.5929 11.85763 2.532708 

dv/dh(xlO A 3) 1.7632 

dr/dh(xlO A 3) 142.1654 

((dh/dt)/dh)*10 A 3)-0.05015 


Altitude 
10 A -3 ft 
********* 
25 

27.63682 

30.43127 

33.36962 

36.43715 

39.61870 

42.89913 

46.26367 

49.69819 

53.18938 


Velocity Radius Onega Tine of Clinb change in 

ft/s ft rad vel Period ft/s altitude 

****************************************************** 
91.44 2970.404 0.030783 221.9309 11.88128 2.636823 
96.08924 3345.269 0.028723 237.8454 11.74902 2.794453 
101.0164 3742.544 0.026991 253.1124 11.60887 2.938350 
106.1973 4160.276 0.025526 267.6375 11.46149 3.067528 
111.6059 4596.372 0.024281 281.3625 11.30764 3.181548 
117.2156 5048.679 0.023217 294.2595 11.14807 3.280427 
122.9997 5515.042 0.022302 306.3255 10.98354 3.364540 
128.9321 5993.363 0.021512 317.5763 10.81479 3.434523 
134.9878 6481.634 0.020826 328.0411 10.64253 3.491190 
141.1435 6977.961 0.020227 337.7583 10.46743 3.535463 


dv/dh(xlO A 3) 3.098 
dr/dh(xlO A 3) 198.2671 
( (dh/dt) /dh)*10 A 3)-0.12209 


Altitude 
10 A -3 ft 
********* 

50 

53.49552 

56.93242 

60.29446 

63.56991 

66.75053 

69.83075 


Velocity Radius Onega Tine of Clinb change in 

ft/s ft rad vel Period ft/s altitude 

****************************************************** 
135.52 6524.540 0.020770 328.9160 10.6274 3.495522 
146.3491 7217.588 0.020276 336.9306 10.20061 3.436899 
156.9966 7899.012 0.019875 343.7327 9.780985 3.362045 
167.4122 8565.595 0.019544 349.5496 9.370495 3.275453 
177.5596 9215.010 0.019268 354.5603 8.970578 3.180611 
187.4131 9845.621 0.019035 358.9067 8.582241 3.080224 
196.9556 10456.32 0.018836 362.7013 8.206161 2.976386 
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72.80714 


206.1765 11046.44 0.018664 366.0344 7.842759 2.870720 


75.67786 


215.0700 11615.61 0.018515 368.9784 7.492258 2.764481 


<3v/dh(XlO''3) 6.72984 
dr/dh(xlO A 3) 681.4859 
( (dh/dt) /dh)*10'3)-0. 29865 


Altitude 
10 A -3 ft 
********* 

75 

77.78992 

80.44966 

82.92280 

85.17866 

87.20602 

89.00718 

90.59302 

91.97938 

93.18450 

94.22732 

95.12639 

95.89920 

96.56185 

97.12890 

97.61334 

98.02662 

98.37880 

98.67862 

98.93366 

99.15047 

99.33466 

99.49108 

99.62385 

99.73652 

99.83210 

99.91316 

99.98190 

100.0401 

100 


Velocity Radius Onega Tine of Clinb change in 

ft/s ft rad vel Period ft/s altitude 

****************************************************** 
212.97 11481.22 0.018549 368.3054 7.575022 2.789921 
231.7457 13382.51 0.017317 394.5158 6.741796 2.659745 
249.6453 15195.09 0.016429 415.8323 5.947448 2.473141 
266.2892 16880.50 0.015774 433.0821 5.208831 2.255851 
281.4707 18417.83 0.015282 447.0373 4.535108 2.027362 
295.1145 19799.45 0.014905 458.3540 3.929625 1.801159 
307.2361 21026.91 0.014611 467.5648 3.391699 1.585839 
317.9085 22107.64 0.014380 475.0931 2.918079 1.386359 
327.2385 23052.43 0.014195 481.2722 2.504035 1.205122 
335.3488 23873.70 0.014046 486.3641 2.144119 1.042822 
342.3668 24584.37 0.013926 490.5756 1.832674 0.899065 
348.4174 25197.07 0.013827 494.0703 1.564163 0.772806 
353.6182 25723.73 0.013746 496.9786 1.333360 0.662651 
358.0778 26175.31 0.013679 499.4052 1.135455 0.567052 
361.8939 26561.75 0.013624 501.4341 0.966102 0.484436 
365.1541 26891.89 0.013578 503.1339 0.821422 0.413285 
367 9355 27173.54 0.013540 504.5602 0.697992 0.352179 
370 ". 3056 27413.54 0.013508 505.7587 0.592812 0.299820 
372.3233 27617.87 0.013481 506.7671 0.503269 0.255040 
374.0397 27791.67 0.013458 507.6162 0.427100 0.216803 
375.4987 27939.42 0.013439 508.3319 0.362350 0.184194 
376.7383 28064.95 0.013423 508.9357 0.307340 0.156416 
377.7910 28171.54 0.013410 509.4452 0.260625 0.132774 
378.6846 28262.03 0.013399 509.8756 0.220971 0.112668 
379.4428 28338.81 0.013389 510.2391 0.187322 0.095579 
380.0860 28403.94 0.013381 510.5464 0.158777 0.081063 
380.6316 28459.19 0.013374 510.8062 0.134567 0.068737 
381.0942 28506.03 0.013368 511.0259 0.114038 0.058276 
381 4864 28545.75 0.013364 511.2118 0.096633 0.049400 
381.216 28518.36 0.013367 511.0837 0.108634 0.055521 
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9.5 V-n Diagram 



Effect of Gust 


100 
110 
120 
130 
1' '■ 
1~_ 
160 
170 
180 
190 
£00 
£10 
£20 
£30 
£40 
£50 
£60 
£70 
£80 
£90 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
4 

41v 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

560 

590 

600 


REM 

REM 

REM 

REM 

REM 

REM 

REM 


REM 

REM 

REM 


REM 

REM 

REM 

REM 


Cent ants 


PI = 3. 141592654# 

PO = (£*PI ) 

Cl = 1.5 
C2 = -1 

11 = . 1 

12 = 10 

Input Varia 


: REM 

Numerical value for 

: REM 

Value for £ * 

pie 

: REM 

max C lift 


: REM 

mir* c lift 


: REM 

Increment for 

Mach # 

: REM 

Increment for 

gust 




INPUT 

"Gross Weight 

" ;W 

INPUT 

"Wing Prea 

” ;S 

INPUT 

"Pspect Ratio 

" ;PR 

INPUT 

"Pltitude of flight 

" ;FP 

INPUT 

"Soeed of sound $ altitude 

" ;SS 

INPUT 

"Density & Pltitude 

" ;R0 


0 |d |d | d I d | d | d l d |d Id I d | d Id lil ( d Calculations I tl Id Id Id | d 0 | d | d | d I d I d I d 

OPEN "es gust si. out" FOR OUTPUT PS 1 LEN=2000 


FOR WG = -40 TO 40 STEP 12 


REM 

REM 

REM 

REM 

REM 

REM 

REM 


PRINT "Wind gust velocity ";WG 

FOR M * . 1 TO 1 STEP 11 
V = M*SS 

P = (00/(1 + (PO/ (PI*PR> > ) > 

N = 1 + ( (RO*P*WG*V)/(£*<W/S> ) ) 

N 1 = (Cl*. 5*R0*(V-''£) / <W/S> ) 

N£ = <C£*. 5*R0*(V''2> / (W/S) ) 

PRINT M ; N1 ; N2 ;N 
PRINT #1, M?N1?N£;N 

NEXT M 


pie 


REM 

REM 

REM 


NEXT WG 
CLOSE 1 
END 
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EFFECT OF WIND OUST 


ALTITUDE SEA LEVEL 


VERTICAL MAGNITUDE 
MACH « 40 20 

******************************** 


0 

0 . 1 
0 . 2 
G . 3 
0 . 4 
0 . 5 
0 . 6 
7 o n o o o 
800000 
900000 
1 


1 

16.91 369 
3 ? . 92733 
48.74107 
64 . 66476 
80 . 56845 
96.48214 
1 1 2 . 3958 
128.3095 
144.2337 
16 0. 1359 


1 

8 . 956844 

16.91 369 
2 1 . 87054 
32 . 82738 
4 0 . 784 22 
48.74107 
56 . 69792 
64.65476 
73.61161 
80.56845 


OF GUST 
1 0 
* 

1 

4 . 9784 22 
8 . 956844 
17.93527 

16.9 1369 
20 . 8921 1 
24 . 87054 
28.84896 
32.82738 
35.80581 
40.78423 


(FT/S) 

-10 


-2 . 9784 2 
-6.95684 
-1 0 .9352 
-14 .9136 
-18.8921 
-22.8705 
-26 . R4 89 
-30.8273 

__ O A p Q 5 O 

-38.7842 


-20 
* : 

1 

6.95684 
14.9136 
22.8705 
30 . 8273 
38 . 784 2 
46.7410 
54.6979 
62.6547 
n 0 . 4 1 1 6 
7 5 . 5 6 8 4 


-4 ; 

v * * * afc * * * . 

-14.913' 
-30.827 
-46.7411 
-62.654 
-78. 568' 
-94 .482 

- 1 1 0 . 3 9 r - 
. - 1 2 6 . 3 0 ( 

- 1 4 1 2 


i n 


O')- 


OUR TO MAX AND MIN CL r S 
MACH n 

************: 


CL = 1.5 


ft : TTTiirr 


5 FA ! 


CL = -1 


0 . 1 

12.0171 

-8.01140 

0 . 2 

43 . 0684 1 

-32 .0456 

0 . 3 

108.1539 

-72 . 1 026 

0 . 4 

192.2736 

-1 28.182 

in 

o 

300.4276 

-200 . 285 

0.6’ 

432 .6 1 57 

-288.410 

0 . 700000 

538.8381 

-392.558 

0.800000 

759.0948 

_ ^ i p t on 

0 . 900000 

073.3856 

-64 8.9 2 3 

1 

1201.711 

-801.140 
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EFFECT OF WIND GUST 


ALTITUDE 25,000 (FEET) 


MACH tt 



VERTICAL MAGNTTIJOF OF GUST 
40 20 10 
************************************ 


(FT/S) 

-10 -20 
*********************'* 


0 . 00000 ! 


1 

7.497005 
13.99401 
20.49102 
26 . 98802 
33.48503 
39 .98203 
4 6 . 

52 
59 
65 


4 7904 
97605 
4 7 3 0 6 
07005 


1 

4 .24 8503 
.7 .497005 
10.74551 
1 3 .99401 
24251 
49 102 
7 3 9 5 2 
98802 
2 3 6 5 3 
4 8 5 0 3 


1 7 , 
20 
23 
26 
3 0 
3 8 


1 

2 . 624 25 1 
4.248503 
5.872754 
7 .497005 
9. 121257 
10.74551 
1 2 .36076 
13.99401 
15.61825 
17.24252 


1 

- 0.62425 
- 2 . 24850 
- 3.87275 
- 5.49700 
- 7.12125 
- 8.74550 
- 10.3697 
- 1 1 .9940 
— 1 3 . 6 ! 9 2 
- 15.2426 


1 

— 2 . 2 4 8 5 0 
- 5.49790 

- 8 . 7 4 5 5 0 
-1 1 . 9940 
- 15.2425 
- 18.4910 
-2 1 .7 395 
- 24 .9880 

_ *■', o O w *- 
-3 1 . 4 850 


-5.4 
-11. 
-18. 
-24 . 
-3 1 . 
-37 . 
-4 4 . 
-50 . 


* * .- 

9 7 ! 

9 9 
4 9 
9 « r 
4 85 
9 8 
4 7 ‘ 
9 7 " 

0 / 



L—O Fact ‘ 22 : CO 

TO MAX AMD 14 1 N 



MACH P 
* * * * * * * * * * * 
0 

0 . 1 
0 . 2 
0 . 3 
0 . 4 
0 . 5 
0 . 6 
n 700 c i) o 

( 1.89 0 0 0 0 
0.900 0 0 0 

1 


►, * * * * * * * * * * * 



CL = -1 

************ * * ****** * * * ; 
n 


4.464943 
17.85977 
40. 1845 
71 .43909 
111 .6236 
160.738 
218.7823 
285 7 : 6 4 

361 .6605 
446.4845 


- 2.97662 
- 11. 9065 
- 26.7896 
- 47.6260 
- 7 4 .4157 
- 107.158 
-145.854 

— 1 n 0 . h 0 '• 
- 24 1.107 
- 297.653 
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tf I- f jr C T 0 E W 1 M 0 GUST 


ALTITUDE 50.000 (FEET) 


ve 


MACH tt 

* * * 3*C * * * 5 
0 

0 . 1 
0 . 2 
0 . 3 
0 . 4 
0 . 5 
0 . 6 
0.700000 
0 0 0 0 0 0 0 


1 


5 

7 


1 3 


PTTCAI MAGNITUDE OF GUST 
40 20 10 

c ****** ********** 

1 1 
.11257 2.056285 
225139 3.11257 

337700 4.168855 
450273 5.226139 

.56285 6.231424 
.67542 7.337709 
7 c< q 8.30300 4 
.00055 0.45 0 2 7 3 


(FT/ 



523 1 4 3 
2.056285 
2 .584427 
3.11257 
3.640712 
4 . 1 63355 

4.606007 
c; o o c i 3 0 


0.47 1857 
-0 .05628 
-0.58442 
-1.11257 

-1 . 64 07 1 
-2 . 1 6835 
_ 2 . 6 0 6 0 0 
_ -> ? 2 5 1 3 


-0.05628 
-1.11257 
-2.76885 
-3.22513 
-4.28142 
-5.33770 
_ ^ 3 6 3 0 
- 7.4602 7 


-1 . 1 1 2 C 

-3.225 

-6.3377 

—7.4502 

— 0 6 6 2 f 
-11.675 

— 13.732 

-15.006 


r\ . j ; ] f i • ') H 0 

1 


! 0 - ? 1 ^ 
.3 1 . 12 5 7 


1 .•? . 5 0 - - - E . 7 5 - ? - "* -'.76. i ' 
11.56 28 5 6 . 2814 26 -4 .28 14'' 


«... O K K r* 

-9 . 5F>?B4 


_ 1 C: n ; ' 

- 2 n . 1 1 c 


LIMIT LOAD FACTOR CONDITIONS 
OIJF TO MAX AND MTM CL'S 


ALT I TUDF 50,000 (FEET) 


M&CH » 


; *• * AT .¥ * 


0 

0 . 1 
0 . 2 
0 . 3 
0 . 4 
0 . 5 
0 . 6 
0.700000 
0.800000 
0.900000 
1 


CL =1.5 CL 

******** ******************** 

0 




1 .383276 
5 . 533105 

> ■}_/[, \ '4 t\ 0 

22. 13242 
34.5319 
49.79796 
67.78055 
88 .52960 
112.0454 
138.3277 


-0 .92213 
- 3.68873 
_ a ; 

- 1 4 .754 0 
_ 1> o _ 0 6 4 6 
-33.1006 

-45. 1870 
- 50.0197 
- 7 4 6 9 6 9 

- 9 2 . 2 1 8 4 
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EFFECT OF WIND GUST 


altitude 


75,000 (FEET) 


VERTICAL 

m in-i tt ^ ® 

************************** 
0 

0 . 1 
0 . 2 
0 . 3 
0 . 4 
0 . 5 
0 . 6 
C . 700000 
O . R00000 
0 . <500000 


1 

1 .64062 1 
? 7 6 12 4? 

2 . 92 1 864 
3.562485 
4 . 203106 
A . 843727 
5 .484349 
6 . 1 2 4 0 7 

6 7655'? | 

7 .406213 


MAGNITUDE 

20 

1 

1 .3203 1 1 
1 ■640621 
1 .960932 
2.281242 
2.601 553 
2.921864 
3.242175 
3.562485 
3.882795 
4 .203106 


OF GUST (FT/S) 

10 _1 0 ~ 20 , 

* * * * * * * * * * * * * * * * * * * * * * * * * ' 


1 

1 .160155 
1 . 3203 1 1 
1 .480466 
1 . 640621 
1 .800777 

1 .960932 

2 . 121087 
2.281242 
2.441308 

7.601653 


1 

0.839844 

0.679689 

0.519534 
0 . 359378 

0.199223 
0 .039068 
-0.12108 
-0.28124 

_0 . 4 4 ’ ? 0 
- 0.60155 


1 

0.679689 

0.369378 

0 . 0 3 0 0 f. 0 

-0 


-0 
-0 
- 1 
- 1 


2 0 . 1 2 4 
6 015 5 
92 186 

? a ? 1 7 
5 6 24 8 
■> 7 * 


0 . 3 
-0 . 
-0 . 
-1 . 
-2 . 

-2 . 
_ O 

-4 . 

- A 


-4 0 

1 

5 9 3 7 f 
? B 1 2 ^ 

g o i a *■ 

5 6 2 4 0 
2 0 3 1 0 
£4372 
4 R4 "« /l 
1 2 '1 0 7 
7 *• 1 

\ 0 5? 


L I M I T 
01 Jr. 


MACH 


LOAD 


AN 


? cor 

M l M 


•• . I T 1 f. 

CL * S 





• * * * * * * ■¥ * "* * * * * * * ^ 
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EFFECT OF WIND GUST 


ALTITUDE 100,000 (FEET) 


MACH tt 

* 

0 

0 . 1 
0 . 2 
0 . 3 
0 . A 
0 . 5 
0 . 6 
0.700000 
0 . 800000 
0 000000 
1 


VERTICAL 

a o 

. * * 

1 

1 .103103 
1 .386380 
1 . 5 7 9 5 8 A 
1 .772778 
1 .965973 
2.159168 
2 .352362 
2.545557 
0 738752 


M A G N I T LI D E 
20 

:*: 

1 

1 . 096597 
1.. 193195 
1 .289792 
1 . 386389 
1 .482987 
1 .579584 
1 . 676181 
1 . 772779 
1 869376 

1 .965973 


OF GUST (FT/S) 
10 


-2 0 —4 C 

:******* *^ 


1 

1 . 04 8299 
1.096597 
1 .1 4 4 8 9 5 
1 . 193195 
1 .24 1493 
1 . 289792 
1 .338091 
1 .386389 
1 .4 34 668 
1 ..482987 


1 

.951701 
.903402 
. 855104 
.806805 
.758505 
.710208 
.661909 
. 6 1 36 1 0 
. 5 6 5 3 1 3 
.517013 


1 

903402 
906805 
7 1 0 2 0 8 

613610 

517013 

420416 

O'llOIS 

C. • 1 ‘ - J 

.227221 

: 3 r i ■") -I 


0.806806 
0.613610 
0.4204 1 c 
0 . 227221 
0.034026 
- 0.15916 
- 0.35236 
- 0.54556 
- 0 . 7 7 f> ■ " 

0 . 9 n 9 ■ 


I i : t i FAC 7 on : ; I ' ■ 

DUE TO MAX AMO MIN CL’S 

MACH It 

************* * * * * * * * **** * * 
0 

0 . 1 
0 . 2 
0 . 3 
0 . 4 
0 . 5 
0 . 6 
0.700000 
0.800000 
0.900000 
1 


CL = 1.5 
* * * - 

0 

0 . 131096 
0.524386 
1.1 7987 
2 .097546 
3.277416 
4.719479 
6 .423736 
8.390186 
10.61883 
13.10967 


I n ) . 0 0 


CL * -1 
********** 

0 

-0.08739 
-0.34959 
-0.78658 
-1 . 39836 
-2.1 84 94 
-3.14632 
-4 .28 2 49 
-5.59345 
-7.07922 
-8 . 73977 
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Appendix A. 10 
Airframe Cost 



The following equations, found in Reference 24, 
allow an estimate to be made for cost associated with 
the airframe only. The equations give an estimate in 
1970 dollars, so the results had to be adjusted to 1990 

dollars. 

.791^,-1.526 . - . 183 

Engineering hours: E= .0396*A *S Q 

where A = airframe weight(lbs) 

S = maximum speed at design 
altitude (knots) 
q = cumulative quantity 
produced 


.873.-1. 89 . _ . 346 

Development support cost: D= .008325*A b y 

.1.160. -1.371 . -1 

Flight test operations cost: F= .01244*A S 

. 764 . 899 178 

Tooling hours: T= 4.0127*A *S Q 

•740 . 543 .524 

Labor cost: L= 28.984*A S Q 


Quality control: Q/C= .13*L 

. 689 . „ . 624*_ . 792 

Material cost: M= 25.672*A S Q 

Solving the above equations for the aircraft, 


.281 


A. 10. 2 



assuming A-3464 lbs (Figure 4.8) and S=346 Knots (Mach 
.44 at 100,000 ft), leads to the following. 

E* ( $ /hr ) +D+F+L+Q/C+M = $10,398,187 (1970 dollars) 


Converting this to 1990 dollars yields a total 
airframe cost of $40,237,700 (assuming a constant 
inflation rate of 7%). 


It is important to realize that this figure is 
probably extremely conservative considering that this 
aircraft is not conventional in any sense. For example, 
the low airspeed (S) and low weight (A) make this 
aircraft appear, to the equations, to be the equivalent 

of a small military aircraft. 
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